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ABSTRACT 

The rocket-borne equipment, experimental method, and data reduction 
techniques used by the Aeronomy Laboratory in the measuronent of electron 
tmperature in the E region are fully described. Electron temperature 
profiles from one daytime equatorial flight and two nighttime midlatitude 
flights are discussed. The last of these three flights, Nike Apache 14.533, 
showed elevated ^--region temperatures which are interpreted as the heating 
effect of a stable auroral red (SAR) arc. 
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i. INTRODUCTION 

The Aeronomy Laboratory is engaged in a wide-ranging investigation of 
the lower ionosphere (50 to 200 km) using sounding rockets. Specific studies 
have included the winter anomaly in the D region, the formation of layers 
(including sporadic-^ layers) by the action of the neutral winds and the 
production and loss of ionization (particularly the effects of energetic 
particles). This report is concerned with one aspect of the program: the 

mr ..suranent and analysis of electron temperature in the E region.* 

Data from three Nike Apache rocket flights are presented. The first, 

Nike Apache 14.532 was launched from Chiles, Peru at 1526 LST (Local Standard 
Time) on 28 May 1975 as part of the NASA launch (^eration Antarqui. The 
second and third (Nike Apaches 14.534 and 14.533, in that sequence) were 
launched from Wallops Island, at 0001 EST on 10 August 1977 and at 0001 EST 
on 5 January 1978, respectively. The value of the 3-hour planetary index of 
geomagnetic activity (Kp) had the value, for the interval including the 
launch, of 2-, 4 and 5, respectively, for the three launches. 

A Nike Apache rocket, of the type used in this investigation, is shown, 
in Figure 1.1, on a military rail launcher at Wallops Island. 

Chapter 2 contains a general discussion of F-region electron temperatures 
during daytime and nighttime and, more specifically, of the role of energetic 
particles and a particular phenomenon: the stable auroral red (SAR) arc. 

The use of Langmuir probes for in situ measurements in the ionosphere 
is described in Chapter 3 in terms of existing theories and regimes of 
operation. This chapter also contains a general description of the probes 
flown on the three rocket flights. 

Details of the implementation of the experirient are given in Chapter 4. 
The equatorial launch (14.532) used both nose-tip and boom-mounted probes 
while the two midlatitude launches used only a nose-tip probe. Chapter 5 
includes procedures for data reduction from digital data tapes. 

The errors in the determination of electron temperature are discussed 
in Chapter 6 together with procedures (mainly in the software) that have been 
developed to minimize them. 

*Although the ionospheric E region is conventionally defined as the altitude 
range of 90 to 160 km we will here, for convenience, use the term to include 
altitudes to 200 km. 




3 


The equatorial launch was planned to allow a eonparison between electron 
temperature measured independently by two rocket>bome probes and by the 
incoherent -scat ter radar at Jicamarca. The results are presented in 
Chapter 7. 

The observations fr<^ the two nighttime midlatitude launches are 
presented in Chapter 8. 

Finally, Chapter 9 contains a summary of the study, some conclusions 
and suggestions for future work. 
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2. ELECTRON TEMPERATURE IN THE £ REGION 
2.1 Heating Sovcposb and Cooling PnooeeaeB 

The ionospheric E region is composed of three interpenetrating fluids, 
these being the neutral gas, the positive ions (atomic and molecular), and 
the free electron gas. Additionally, in the lower E region, there are small 
numbers of negative ions formed through the process of electron attachment 
to neutral species. These negative ions are of no iiq)ortance in terms of 
the ^-region thermal balance because of their low concentration and will be 
ignored. Each of these three fluids has a characteristic ten^erature, a 
specific heat, and is in thermal contact (tc a greater or lesser degree) 
with the other two con^nents through collisional processes. At all 
altitudes and under all conditions most of the thermal energy is contained 
in the neutral gas. 

2.1.1 Lonnef E region. At altitudes less than about 130 km there are 
sufficient collisions that the three fluids have the same temperature 

. Evidence for this position has been obtained by incoherent - 
scatter radar studies of the lower E region by many researchers, for example: 
Carru etal., 1967a, b; EvanSt 1967; Wand and PerkinSt 1968, 1970; and 
Salah et at., 1975. At these altitudes the electron temperature will 
immediately follow any changes in the neutral gas temperature, and the 
question of heating and cooling the electron gas reduces to asking what 
processes heat and cool the neutral gas. 

During the day solar ultraviolet radiation in the range 1026-1750 X 
is absorbed by molecular oxygen causing photodissociation. The atomic 
oxygen produced in this manner liberates energy mainly as heat upon 
recombination. Solar X-rays (X > 1 nm) are also important. 

Another source of heating between 95 and 140 km is the dissipation of 
internal gravity waves, with periods between 5 minutes and 10 hours, 
propagating upward from the mesosphere. These waves yield energy to the 
neutral atmosphere up to heights of about 110 km by crijating turbulence 
which eventually degrades to heat. At higher levels turbulence ceases 
[Blamont and de Jager, 1961]; the upper boundary of this region of mixing 
is called the turbopause. 

Above the turbopause the dissipation of gravity-wave energy is due to 
molecular kinematic viscosity. The estimated heating rate due to these 
waves at midlatitude are 10 K/day at 95 km, 30 K/day at 105 to 110 km 
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and 100 K/day near 140 km [Uineet 196S] . If Hines is correct in these 
estimates, then the heating due to gravity-wave dissipation is cranparable 
to heating due to solar-radiation absorption in the lower E region. Hines 
has put this forth as an explanation of the rapid rise of temperature (with 
increasing altitude) that occurrs in the E region. 

Another major heat source is adiabatic compression of the atmosphere 
by tidal motions. The tides are thought to be driven by the heating of 
stratospheric-mesospheric ozone by UV radiation during the day [Butler and 
Small, 1963; Lindzen and Chapman, 1969] and are often called thermal tides 
to emphasize this point. Vtost of the energy of the ozone driving mechanism 
goes into the first semidiurnal (2,2) tidal mode causing an upward and 
downward propagating 12-hour pressure oscillation. The downward propagating 
wave continues to the earth's surface and is observed there as a semidiurnal 
pressure oscillation. The upward traveling (2,2) wave does not propagate to 
the E region, becoming evanescent below the mesopause (at about 80 km) . 
However, there is mode mixing in the mesosphere, coupling the (2,2) mode 
energy to the (2,3), (2,4) and (2,5) semidiurnal tidal modes. The (2,4) 
mode, having a smaller equivalent depth than the (2,2) or (2,3) modes, can 
propagate into the lower E region. As the wave passes through a given 
altitude level, it causes a 12-hour rise and fall of the neutral gas 
temperature T due to adiabatic compression and expansion. This modulation 
of the temperature has been observed by the incoherent-scatter radar 
facilities at Millstone Hill (42* N) and at St. Santin (45* N). It typically 
amounts to an 11% variation in the mean temperature at 116 km but on occasion 
as much as 20% variation is found [Salah and Wand, 1974]. 

There are other sources of heat in the midlatitude lower E region which 
usually are not comparable in importance with the sources already discussed. 
These include forms of chemical heating (other than the oxygen heating), 
joule heating due to ionospheric dynamo electric currents, and heating due 
to energetic-particle precipitation during times of moderate to severe 
magnetic disturbance. This last effect might very well be more important 
than previously assessed, particularly as an ionization and heating sources 
at night. The Aeronomy Laboratory is currently engaged in a study of 
energetic-particle flux in the E region using rocket-bome solid-state 
particle detectors and data available from several orbiting satellites 
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[AtitTi «t al,, 1974; Voas and Smith, 1979]. This source will be discussed 
separately in Section 2.2. 

In all cases discussed here the neutral gas is colled by downward 
conduction of heat to the mesosphere. Above the turbopause this heat 
conduction is acco^lished by molecular diffusion. At the turbopause and 
below the mechanism for heat conduction is eddy diffusion which carries 
heat through the temperature minimum at the mesopause and deposits the 
energy finally at heights of about 50 km. 

Other sources of cooling such as radiative cooling by photoemission and 
chamical cooling by transport of atomic oxygen exit of the region of 
production may also be significant. 

2.1.2 Upper E region. At these altitudes the principal source of 
energy is also solar extreme ultraviolet radiation (EUV) . At midlatitudes 
during the day there are no other competing sources of heat. Photoionization 
occurs, peaking at about 150 km, mainly due to absorption of radiation in 
the spectral range 140-796 X [Sohunk and Walker, 1973]. The energy of the 
ionizing photons well exceeds that required for ionization, the excess 
energy .oing into ion and electron kinetic energy, and electronic excitation 
of the resultant ion. The ion electronic excitation energy is typically less 
than 4 eV, leaving some tens of electron volts to be transferred to ion and 
electron kinetic energy. Most of the kinetic energy is carried by the 
electron, because of the large mas ratio between the ion and electron, 
resulting in hot photoelectrons, '"^ese photoelectrons provide the entire 
heat input to the ambient electron gas. The route by which the photoelectrons 
share their energy with the population of thermal electrons is complicated, 
the general outline being given in Figure 2.1. Each photoelectron gives up 
most of its energy by effecting impact ionization of ambient neutrals, 
creating other excited ions and hot electrons. This will happen repeatedly 
and may be thought of as an energy cascade, since each impact ionization 
will remove 25 eV (mean) [Pharo et al., 1971] from the photoelectron. 

Between impact ionizations the photoelectron loses energy in a continuous 
manner by elastic Coulomb collisions in the ambient thermal electron gas. 

The many hot electrons due to impact ionization also yield energy to the 
thermal population by Coulomb scattering, which is highly efficient because 
the masses of the scattering particles are equal. At altitudes less than 
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200 km heating the thermal electron gas in the above manner is strictly a 
local heating process, meaning that the energy of the photoelectron will be 
dissipated within one neutral scale height of its origin site. 

Solar EUV also causes direct heating of the neutral gas by formation 
of excited atcxBS and molecules. This energy is degraded and shared with 
the ambient neutral and ion gas by st^er-elastic collisions. Another less 
inqiortant source of heat for the neutral gas is excited photoions and hot 
ions created by photoelectron impact. If a general accounting of the heat 
input to the atmosphere is done, it is found that most of the heat input 
from solar EUV goes directly into warming the neutral atmosphere - the heat 
input to the neutral gas is roughly an order of magnitude greater than the 
heat input to the electron gas [Sohmk and Walker ^ 1973) . However, the 
specific heat capacity of the electron gas is several orders of magnitude 
smaller than that of the neutral gas. It is expected, then, that the 
electron temperature T should be greater than the neutral gas temperature 

T , provided that the thermal contact between the two gases is weak enough 

yi 

to permit this to occur. This is precisely what happens, as satellite 
measurements in Figure 2.2 demonstrate. The lack of thermal equilibrium 
between the electron and neutral gas is the most striking difference between 
the ionosphere above 130 km and that below 130 km. Thermal nonequilibrium 
was first indicated by the rocket measurements performed in 1958 by Boggeee 
et al. [1959] and verified since then by rocket, satellite, and radar 
studies. Energy flow from the hot electron gas to the warm neutral gas is 
moderated by the rather poor thermal contact between the two gases due to 
the neutral 'electron collisional frequency decreasing with height. 

Positive ions receive kinetic energy during the day in the ionization 
process and from Coulomb interactions with hot electrons. However, the ions 
are immersed in a sea of neutrals for which energy transfer by collision 
is a most efficient process, again, because of equal masses. This excellent 
thermal contact constrains the ion temperatures to be approximately equal 

to each other and equal to T . Often the ion temperatures for specific 

■f ” 

species are all taken to be the same and are simply given as T^. During 

the day (T. •• T ) is usually no more than 10 K at 200 km, although sometimes 

V ft 

it is as much as 100 K. A recent effort to measure T , T., and T 

fl % e 

simultaneously from the Atmosphere Explorer-C (Explorer 51) satellite is 
described in Brace et al, [1976]. Figure 2.2 is from this work, showing 
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T , T., T and the ion concentration near perigee for a midlatitude pass 
during February. At 200 la the latitude was 33* N, the local time 1200 hr, 
and the solar zenith angle 48*. The pass occurred under moderately quiet 
geomagnetic conditions, the three-hour geonagnetic planetary index (Kp) 
being 24. At 160 km the ion and neutral gas temperatures are the same at 
approximately 700 K. The electron temperature at the same altitude is 
about 1150 K. It is also clear from the diagram that in the F region the 
full heating and cooling theory for ions must be considered, as there is a 
significant departure from thermal equilibrium between ions and neutrals. 
However this is above the region of interest of this report. 

The electron gas is cooled by a variety of processes which are indicated 
in Figure 2.1. Excitation of the fine structure levels of atomic oxygen is 
the dominant cooling mechanism for the electron gas. The fine structure 
levels are so closely spaced that they are easily excited even by thermal 
electrons, the energy then going into infrared radiation at 63 ym. Also 
important is the rotational excitation of N 2 and O 2 . At elevated electron 
temperatures, say greater than 1500 K, cooling of the electron gas due to 
vibrational excitation of N 2 and O 2 must also be considered. Table 2.1, 
from Sohutz et at. [1975], is a list of quantitative cooling formulas 
coimonly used in modeling the thermosphere. 

The ultimate fate of the electron gas thermal energy is to be 
transferred to the neutral gas, or to be lost as radiation at wavelengths 
for which the thermosphere is transparent. The former process has a small 
but significant effect on the neutral atmospheric temperature. Heat input 
to the neutral gas from the electron gas comprises roughly 15 percent of the 
total neutral gas heat influx {Hertnan and Chandm, 1969a, b]. Once energy 
is absorbed by neutral atmosphere, it is conducted downward by molecular 
diffusion or is lost to atomic oxygen fine structure excitation with 
subsequent emission at 63 ym. 

2.2 Energetio Fartiale Precipitation and SAP Area 

The two preceding sections addressed the thennodynamics of the 
midlatitude E region during daylight hours. Conditions at night are 
markedly different. Without the ionizing radiation from the sun the electron 
concentration drops by roughly a factor of 100 (on a magnetically quiet 
night) within two hours after sunset. In the absence of solar radiation the 
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T«bl« 2.1 Quantimive electron cooling foraulee used In aodellng the 
thermosphere. Cooling rates are in units of (eV cn'^s*^). 
Cross section calculations by Ho«gy [1976] suggest the atomic 
oxygen fine structure cooling coefficient should be reducsd to 
roughly 60% of that shown. 


Atomic oxygen fine structure excitation [Dalgcamo, 1969] 


- 12 , 


I«][0](3.4 K 10 “)(!T^ - yci.O 


7.0 X 10’® T^) 


n 


N 2 rotational excitation [Dalgaxmo, 1969] 
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O 2 rotational excitation [Dalgamo, 1969] 
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N 2 vibrational excitation [Dalgamo, 1969] 
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solar theroal input to tho electron gas is zero, and the electron tcaporature 

falls to the tenperature of the neutral gas. Representative values for an 

5-3 3 

altitude of ISO kn are n equal to 3 x 10 cn during the day and 3 x lo 

-3 * 

cm at night, and equal to 700 K during the day and 500 K at night. 

The nighttime electron concentration and temperature are strongly 

dependent on magnetic activity; presently this is thought to be due to 

energetic -particle precipitation. Nike Apache 14.439 was launched from 

Wallops Island on the night of a severe magnetic storm (Kp • 8; 0003 EST 

on 1 Novanber 1972) and carried a Geiger counter in the payload. Energetic 

-2 -1 

particles (>70 keV) \ere detected with a maxinaiffl flux of 1086 ± 26 cm s 

ster'^ at 140 km . et aZ., 1974]. The energetic particles were 

presumed to be electrons until subsequent experiments were carried out using 
solid-state surface-barrier detectors with broom magnets. Instrumentation 
such as thi:; aboard Nike Apache 14.533 (Kp ■ 5; 0001 EST on 5 January 1978) 
hris identified the flux of particles as consisting of protons and possibly 
0* ami lie* Ions \Vo&b and SHth, 1979]. 

Usinf^ electron -concentration profiles from nine nighttime miulatitude 
launches, Voae and [1979] have calculated the corresponding ionization 

?ates. These rates aru plotted in Figure 2.3 as a function of magnetic 
activity. The solid line is a least-squares fit, with a correlation 
-roefficient of 0.93 for the nine observations. The ionization rate appears 
to be causally related to the level of geomagnetic activity. Present data 
indicate that the upper limit for all ionization sources (other than 
energetic-particle fluxj in the upper E region is lO”^ cm"^ s“^. This 
supports the conclusion that energetic-particle precipitation is the 
principal source of ionization in the midlatitude E region on geomagnetical ly 
ouiet rights as well as disturbed nights. (Even under the most disturbed 
conditions particle precipitation is insignificant as an ionization source 
during the daylight hours, amounting to less than 1 percent of the daytime 
ionization rate.) 

With the nighttime ionization rate controlled by geomagnetic activity, 
it is clear the electron tem>erature will depend on Kp, for each particle 
ionization event ultimately puts about 35 eV of energy into the f-region 
electron gas. Bnith et al. [1974] have analyzed electron tenq)erature data 
from Evans [1973] at Millstone Hill and found the electron temperature at 
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Figure 2,3 The variation with laagnetic activity of the 
innjjation rate in the upper E region near 
midnight [Voea and &nith, 19791. 
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200 kn 1* iftll represented by (112 ICp * 588} X for e renge in Kp of !• to 
S>. The correUtion coefficient is O.M. 

Nighttino particle precipitation can qvialitatively explain midlatitude 
F>region electron tenperatures on quient and aoderately disturbed nights. 
However , during geomagnetic storms with a Kp of 5 or greater another heating 
nwchanisn s^etimes takes over, this being the stable auvoral red (SAR) arc. 
The SAR arc (or M-arc, for midlatitude arc) was first observed in the night 
sky from southern France by Barbier in 1956 [Barbier, 19SS] . The arc is a 
subvisual red aurora of amazing spectral purity; the wavelength of 6300 X is 
the emission due to the atomic oxygen translvion 0 I(^p 2 - ^^2^* absence 
of other spectral lines, most notably the 5577 I line of atmnic oxygen (sec 
Figure 2.4), implies low energy excitation of the oxygen atom since only 
1.97 eV are required for the ^0 transition. 

The arcs run east to west and extend around the night>sidc of the earth. 
The arc of September 1967 was viewed by the QGO 4 satellite and shown to be 
globe encircling [Feed and Blmont, 1968]. Generally, SAR arcs extend 
several hundred kilometers in the north-south direction. The position of 
SAR arcs is usually between the magnetic L shells 2 and 4 [Roach and Roach, 
1963] . 

The mean altitude of the most intense 6300 X emission has been found 
to lie near 40C km [Roach and Roach, 1963]. Tohnatsu and Roach (1962] have 
produced a normalized isophot*.^ plot of a typical SAR arc in Figure 2.5. 

The occurrence of SAR arcs more than an occasional phonomenon. Their 
frequency follows the solar cycle; during the last solar cycle period SAR 
arcs were reported on 20 percent of the viewing nights at Fritz Peak, 

Colorado [Marovtoh, 1966]. 

There are two major atmospheric properties associated with SAR arcs: 
the arcs form within or on the equatorward edge of a general ionospheric 
electron-concentration depression. This is shown dramatically in Figure 2.6 
[Clark et at., 1969] where the location of the arc, as determined by ground- 
based airglow observations, is indicated by the blcck bar. The electron- 
concentration data were obtained from the topside ionospheric sounder aboard 
the Alouette II satellite. 

Secondly, there is a strong enhancement of the electron temperature 
CMipared with regions outside the arc [Horton and Findlau, 1969; Chayidm 
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Figure 2.5 Normalized i 6300 X isophotal representation 
of an SAR arc cross section [To'lvnatsu and 
Roach, 1962]. 
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et aU, 1971; nolle et al. , 1971; Nagy et al,^ 1972], Satellite temperature 
measuranents above SAR arcs are shown in Figure 2.7 [Roble et al.^ 1971]. 

In each case the electron temperature peaks just above the photometric 
location of the arc. These measurements were taken at roughly 1000 km, but 
a marked warming would also be expected in the upper E region due to downward 
heat conduction. One of the launches analyzed in this report (Nike Apache 
14.533) is believed to have flown beneath an SAR arc and detected elevated 
electron temperatures due to this mechanism. 

SAR arcs occur on magnetic sheets which pass through the plasmapause 
region [Carpenter, 1971; Chappell et al,, 1971; Hooh and Smith, 1971]. 

Given the worldwide character of SAR arcs, most researchers agree that 
heating does not originate in or immediately above the ionosphere, but takes 
place in the magnetosphere with the proton ring current as an energy source 
as shown in Figure 2.8 [Reee and Roble, 1975]. This presents a difficult 
transport problem: how does the energy reach the ionosphere in sufficient 

amounts to excite the arc and cause general heating, while not providing too 
much energy (which would excite other spectral lines than just the oxygen 
6300 A line)? Presently there are four hypotheses which separately or 
collectively may explain SM arcs. All agree on the necessity of a 
geomagnetic disturbance initally. And all agree that ultimately there are 
electron energies of 2 cV in the upper F region, with consequent collisional 
excitation of atomic oxygen to the state, and spontaneous emission at 
6300 X. They differ regarding the transport mechanism. 

Cole [1975] claims Coulomb interaction between plasmaspheric electrons 
and energetic ring particles (H*, He*, 0*) is the dominant heat source. The 
electrons then carry their heat to the ionosphere by traveling along magnetic 
flux tubes due to the high thermal conductivity along field lines. 

Reee and Roble [1975] believe that intense ion-cyclotron waves are 
generated in the equatorial region due to the overlap between the ring 
current and the plasmasphere. These plasma waves give up their energy to 
plasmaspheric electrons by Landau damping. These hot electrons then stream 
down into the ionosphere along field lines, as in Cole's hypothesis. 

The hypothesis of Cornwall et al. [1971] is that the initial interaction 
in the equatorial ring current in Coulomb interaction, warming the ambient 
electrons. After some warming ion-cyclotron wave generation begins. The 
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Figure 2.7 The electron temperature and electron density :.ieasured by 
the Alouette 2 satellite as a function of the magnetic 
invariant latitude. Curves A-E were measured at different 
times during the geomagnetic storm of October 29-November 2, 
1968. The hatched areas indicate the calculated X65C0 
emission rate in Rayleighs, and the solid areas indicate 
the X6300 emission rate measured from airglow observatories 
\Eohle et al. , 1971] . 
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energy from the Ion-cyclotron waves is lost to the electrons by Landau 
damping again. Cornwall et al. argue that Landau damping is not effective 
until the electron gas has been preheateu by the Coulomb interaction. When 
Landau damping becomes effective they call it Landau resonant energy 
exchange. Once the electrons are heated they carry the energy to the 
ionosphere along magnetic field lines. 

HasegcoM and Mima [1978] have an hypothesis radically different from 
the previous three. They argue that due to two processes 
(magnetohydromagnet ic surface wave generation and plasma drive wave 
instabilities) kinetic Alfven waves are created. This is a type of shear 
Alfven wave with an electric field parallel to the geomagnetic field. 
Electrons are heated by resonant energy exchange at the plasmapause and 
transported to the ionosphere by the parallel electric field. They show the 
kinetic Alfven wave wavelength is comparable to the length of a geomagnetic 
field line and that such long distance transport of hot electrons is quite 
reasonable due to the parallel electric field. 
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3. LANGMUIR PROBES ON SOUNDING ROCKETS 
3.1 Plama in the Region 100 to 200 hn 

The region 100 to 200 km has electron concentration of the order of 
10^ cm*^ during the day. The neutral particle concentration at 150 km is 
about 5 X lO^^cm"^ (see Table 3.1) so that the region is only slightly 
ionized; a very tenuous plasma. In the context of plasmas generated in 
the laboratory the ionosphere may be considered quiescent since the 
characteristic time of change is of the order of an hour if we exclude the 
transient phenomena occurring at sunrise, sunset, and during solar eclipses. 

Because of the hypersonic velocity of sounding rockets, typically 
1.5 km s~^ at 100 km altitude, the positive ions and neutral gas impinge 
upon the rocket with a well directed flow velocity. However, the electron 
mean thermal velocities are 200 to 300 times greater than the ion thermal 
velocities, meaning that it is safe to assume the main electron population 
has an isotropic velocity distribution relative to the rocket vehicle. As 
Boyd [1968] has noted "... the supersonic or hypersonic velocity of the 
spacecraft gives rise to a situation in which the ion mean energy exceeds, 
by perhaps over an order of magnitude, the mean energy of the 
electrons — a complete reversal of the situation to which the classical 
Langmuir papers are addressed." A further complication resulting from the 
high rocket velocity is local aerodynamic heating, though this is limited 
to the lower E region. 

A plasma is described by several characteristic distances, the most 
important of these being the electron Debye length This length, also 

called the Debye shielding length or (erroneously) the plasma sheath 
thickness, is the approximate range of Coulomb forces in plasma. Normally, 
one thinks of the Coulomb (electric) field due to a net charge as having an 
infinite range. This is not the case in a plasma where the free electrons 
(and to some extent, the ions) redistribute themselves around a net charge 
so as to shield the electric field from the distant plasma. Thus, when a 
conducting body is placed in a plasma it will generally have a space 
charge layer, the plasma sheath, surrounding it. Within the sheath charge 
neutrality does not hold, allowing electric fields to be maintained across 
the sheath. 


Table 3.1 Some nominal daytime values of basic plasma parameters in near space 
for mean solar activity [Boyd, 1968], Research in the past decade 
indicates T = T - T at 100 km. 
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The Debye length may be considered in a rough sense as the electrical 
"range of disturbance" in plasma due to a rocket body or electric probe, 
although it is important to keep in mind that the Coulomb field does not 
suddenly drop to zero at a distanct of X^. By placing a sensing probe more 
than several Debye lengths from the rocket body, for instance, one can be 
reasonably sure that the probe is immersed in the ambient plasma. This is 
with the provision that the neutral gas dynamics does not cause the probe to 
be in a region of rarefaction (the wake). 

The electron Debye length is given by the formula 

X^ - (3.1) 

in c.g.s. units, where k is Boltzmann's constant, T is the electron 

® -3 

temperature in Kelvin degrees, is the electron concentration in cm , 
and e is the electronic charge. Debye lengths for different altitudes of 
interest are given in Table 3.1, in «;>ach case being less than 1 cm. 

Following the Debye length in inoortance is the neutral-neutral or 

ion-neutral mean free path X^ (th«> are assumed equal), and the electron-ion 

mean free path X^^. These lengths are indicative of the collisional coupling 

between the plasma particles and are important in the following sense: if 

the mean free paths are much greater than the plasma sheath thickness, then 

any disturbance to the plasma due to particles emitted by the probe surface, 

reflected from the probe surface, or absorbed by the probe surface will be 

slight due to the large volume of plasma involved. Large volume here means 

that the collisional volume is much greater than the volume of the plasma 

sheath. Since the response of an electric probe primarily depends o'l the 

plasma within the sheath, the measurement errors arc minimized if X and X 
r t o e* 

are much greater than X^. This situation is usually referred to as the thin 
sheath case. 

From Table 3.1 we see that X and X are much greater (by factors 

o e* 

>20,000) than X^ at 150 and 200 km. However at 100 km X^/X^„ - 100. Although 
this is still considered to be a thin sheath case, some caution in the use of 
conventional Langmuir probe theory is warranted, for to an increasing extent 
collisions are occurring within the plasma sheath. As alluded to earlier, 
there is the possibility in this situation of collisional heating within 
the sheath volume. 
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Before leaving this section it is necessary to emphasize that the most 
important plasma parameters are the ion and electron temperatures and the 
electron concentration. The Debye length and mean free paths have been 
discussed first because they are important concepts in understanding 
qualitatively how errors are introduced in probing plasma with a practical 
probe. 

The electron temperature is the quantity which (to those not working 
with plasma) is most alien to our senses. For example it is initally 
surprising to find that the electron temperature within a typical fluorescent 
lamp is well in excess of 20,000 K. First, then, we must separate our 
association of heat with temperature, for it requires very little energy 
to greate very high electron temperatures in a plasma. This is due to the 
very small heat capacity per unit volume of electron gas, which is given by 

(Z/2)kn (c.g.s. units). The heat capacity of one cubic kilometer of the 

® 5-3 4 

£■ region during the day (with n « 10 cm ) is 2 w lo ergs/K; which is 

® 4 

comparable with the heat capacity of the head of a pin (3.6 x 10 ergs/K 
for one cubic millimeter of iron). 

Energy deposited in the midlatitude E region by energetic-particle 
precipitation does not cause any perceivable effect during the day. 

However, when the heat capacity decreases at night the particle energy flux, 
although small, is reflected in an elevated electron-temperature profile. 
Clearly, the electron temperature is a very sensitive ionospheric parameter 
and can be affected by even minor sources of energy. 

The mean kinetic energy of an electron in the population is given by 

~m = (1.293 x 10“'*) T (cV) ':.3.2) 

2 e e 2 e e 

Sometimes is called the electron kinetic temperature, emphasizing this 
point and distinguishing it from other possible definitions of temperature. 

A 1000 K plasma, then, has a mean electron energy of only 0.13 eV and an 
r.m.s. electron velocity of 213 km s If a Maxwellian velocity 
distribution is appropriate, as it is in the case of electronic thermal 
equilibrium, then the electron temperature conveys much more information. 
Figure 3.1 gives the velocity distribution for a 1000 K Maxwellian plasma. 


ENERGY (cV) 



tribution for a 1000 K Maxwellian plasma 
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3.2 Langmuir Probta 

The Langmuir probe (or electrostatic probe) was among the first 
experiments to be carried into the upper atmosphere when the American rocket 
program got underway in 1946 using Gexman V>2 rockets. The probe is simple, 
consisting of a conducting material that is inserted into the plasma on an 
insulating mechanical support. The potential of the probe is varied with 
respect to the rocket body and the resulting current flow to the probe is 
recorded on board the vehicle or telemetered to ground on a telemetry link 
of sufficient bandwidth. The resultant current/voltage curve, the prone I-V 
characteristic, depends mainly on the ion and electron temperatures and the 
plasma concentration and, to some extent, on the magnetic field, the flow 
velocity, and the condition of the probe surface. 

The Langmuir probe is capable of making direct local measurement of 
plasma parameters. This is of great importance in the study of acronomic 
phentxnena which may be transient or small-scale, for other rocket and ground- 
based techniques involved spatial or time integrating (averaging) methods 
necessitated by signal-to-noise considerations. If care is taken in the 
design of an electrostatic probe and its associated electronics package, 
spatial and time resolution of 1 m and I ms can be achieved from a vehicle 
traveling 1 km/s. 

The essential simplicity of the Langmuir probe is scxnewhat offset by 
difficult theory and a variety of experimental problems which may beset the 
experimenter. Probe theory is complicated because the probe surface is a 
boundary to the ambient plasma and near boundaries the equations that govern 
plasma behavior change. There is no general theory for the Langmuir probe 
response under arbitary plasma conditions, and many of the experimental 
problems have yet to be adequately addressed. However, the promise of the 
Langmuir probe (direct measurement of local properties) is great and the 
theory is relatively simple over the altitude range 100 km to 800 km. 

3.2.1 Probe-operatiorial regvnea. Probe-operational regimes are usually 
described with the ratios ^nd where /? is the characteristic 

dimension of the probe, ^ , is a representative mean free path, and ^ r, is the 
Debye shielding length. These ratios determine the domains of probe 
operation, which are: collisionless/collisional , thin/thick sheath, orbital 

motion limited current; etc. The various probe-operational regimes arc 
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shown in Figure 3.2. along with a curve indicating the various plasma 
conditions which are attained in the ionosphere. 

Most of the D region of the ionosphere (specifically the range SO to 
8S km) is contained in the collisional-sheath domain. Here the current to 
the probe is controlled by collisional processes occurring within the plasma 
sheath. This is the problem of a flowing continuum low-density plasma which 
Chcmg ccnd Laftvffnboiia [1976] have numerically solved for specified plasma 
conditions and probe geometries (sphere in flow, and cylinder in cross flow). 
As might be expected the solutions in this region depend strongly on probe 
geometry, vehicle velocity, and (.for the cylindrical probe) angle of attack. 
For given plasma conditions and flow velocity Chang and Laframboise were 
able to calculate the total probe current, each calculation for a given probe 
potential requiring 20 minutes of numerical computation on a CDC 6600 
cominiter. While it may be argued that it should be possible to solve the 
reciprocal problon (that of determining 0-region plasma parameters from a 
Langmuir probe characteristic using first principles), it is not reasonable 
to do so at the present. Furthermore, a sensitivity analysis might indicate 
that this approach would never be worthwhile as a research technique. For 
example, at an altitude of 50 km the Debye length (A^) might be 40 cm, 
meaning that the plasma sheath due to the rocket body could be over a meter 
in thickness. It is a simple enough matter to create an isolated probe in 
a digital computer model, but quite a different affair to do so in the 
D region from a small sounding rocket. 

It should be noted that the present use of Langmuir probes in the 
r region is restricted to a technique developed by Smitn [1963] for measuring 
fine structure (on a scale of meters) in electron concentration. No attempt 
is made in this method to determine electron concentration from probe 
current using first principles. Rather, the probe current is recorded as a 
function of altitude with a constant positive potential on the probe. The 
assumption is made that the probe current and electron concentration arc 
related by a constant factor, n /t which changes only slowly with altitude. 

'S 

This calibration factor is obtained fron an independent experiment, 'isually 
a radio -propagation experiment. 

Referring to Figure 3.2 again, at about 85 km is the transitional sheath 
probe regime. Some effort has been made to deduce electron temperature from 
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probe measiironent in this region, although the results have not been 
encouraging. Ion-neutral collisions are still niunerous enough to render 
the conventional Langmuir probe theory invalid. 

Extending from approximately 100 km to in excess of 800 km is the region 
where conventional Langmuir probe theory is valid. This is the thin sheath 
region, mentioned in the previous section. Here the theory of Langmuir probe 
operation is based on the kinetic theory of gases with the additional 
assumption of there being no collisions within the plasma sheath. This is 
an excellent assumption and results in the elegant and powerful theory, 
described in detail in the following section (collisionless probe theory), 
which is appropriate to most of the work described in this report. 

Beyond 800 km the Debye shielding length begins to increase due to the 

high electron temperature found in this region and the gradual decrease in 

electron concentration (above the F-region peak). As increases, and in 

the complete absence of collisions, particles which are attracted to the 

probe tend to execute orbital motions about the probe for particular values 

of energy and angular momentum. This increases the residence time of the 

attracted particle near the probe. However, as in the thin sheath regime 

at lower altitudes, the electron temperatare is determined by analysis of 

that portion of the probe characteristic where electrons are repelled from 

the probe surface (retarding potentials). The retarding potential analysis 

developed in the following section is applicable for the determination 

of T in the orbital-motion-limited case also. 
e 

3.2.2 Collisionlese Langmuir probe theory. The plasma sheath, in 
addition to being dependent on the various plasma parameters, is a function 
of the probe potential. There is a potential for which ions and electrons 
adjacent to the probe are neither attracted to nor repelled from the probe 
surface. This potential is called the plasma potential or space potential. 
There is no space charge sheath surrounding the pi'obe at plasma potential. 

If one were to microscopically examine the plasma near such a probe, one 
woudl find the distribution of speeds (for both positive and negative 
particles) is identical with that in the ambient plasma. Adjacent to the 
probe surface the particles are directed only towards the probe surface 
(assuming the surface is perfectly absorbing and nonemitting) whereas the 
particle velocities are isotropic in the ambient plasma, but this is of 
minor importance. 


31 


When the probe is at plasma potential, the ion and electron currents 
to it are controlled strictly by the thermal velocities of the ions and 
electrons in the surrounding plasma. From the kinetic theory of gases the 
number of particles striking a unit area per second is ny/4 where n is the 
particle concentration and u is the mean particle velocity. The magnitude 
of the electron and ion current density to the probe surface is given by 





(3.3) 

and 

U'.l 

= n ey /4 
+ + 

(3.4) 


The assumption that the probe is perfectly absorbing (fully 
accommodating) is not always satisfactory, in which case it is necessary to 
reduce the current densities in equation (3.3) and (3.4) by experimentally 
determined accommodation factors ranging between 0 and 1. The assumption 
that the probe is nonemitting is very good. Probes can be made to emit by 
applying excessive potentials or by heating the probe (thermionic emission) , 
but that is not the intent of the present work. 

The electron temperature may exceed the ion temperature because 
the energy of the photoelectrons is given to the electron gas. For the 
moment, however, assume that the ion and electron gases are in thermodynamic 
equilibrium. The ion mean energy is then equal to the electron mean energy 

hm V ^ = Hm (3.5) 

e e * * 

The ratio of electron to ion mean velocity is, therefore, 

- - - l' 

u /y = (m /m ) ^ (5.6) 

e + ^ + e 

or about 170, depending on the ion composition at a given altitude. 
Throughout the E region there is strict charge neutrality, so that, = n^. 
This fact, together with equation (3.6) above, indicates that when a probe 
is at plasma potential the electron current to its surface will exceed the 
ion current by a factor of about 170 (greater in reality, where there is 
not thermodynamic equilibrium). Clearly, the probe must take in a 
considerable electron current if it is to remain at plasma potential. 

If measurements are to be made on a rocket with an electrostatic probe 
at plasma potential, it will be necessary to provide a means to balance the 
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net current which the probe collects from the plasma with an equal but 
opposite current elsewhere in the system. As mentioned earlier, plasma 
potential is a convenient potential to use since the current densities 
are immediately obtained from kinetic theory. Therefore, the return current 
is an important consideration. 

Another special case occurs when the probe is biased so as not to draw 
any current. Since the ion and electron currents are then equal in magnitude 
(but of opposite sign) it follows that the ion and electron velocities 
adjacent to the probe are perturbed relative to the ambient plasma. This 
point may be seen by considering a probe, initially at plasma potential, 
which is then isolated electricc) ;>• irom the spacecraft. The electron 
flux to the probe initially exceeds the ion flux, and the probe quickly takes 
on a negative potential of such a magnitude as to make the ion and electron 
currents equal. In this situation the probe is surrounded by a negative 
plasma sheath which acts to accelerate positive ions to the probe surface 
while retarding electrons. The probe is at floating potential. It is 
relative simple experimentally to determine the probe floating potential 
with high resolution by sweeping the probe voltage with respect to the space 
vehicle and recording the applied voltage at which the probe current is 
precisely zero. The plasma potential is more difficult to measure. 
Fortunately it is possible to determine the electron temperature without an 
accurate knowledge of the plasma potential. 

Since the probe electron current is more than two orders of magnitude 
greater than the ion current when the probe is biased near plasma potential, 
it is worthwhile initially to consider all the current as electron current. 

and Langmuir [1926] developed the theory for the current to a 
probe in a Maxwellian plasma. For a probe of any geometry as its voltage is 
made negative with respect to plasma potential the electron current density 
will decrease exponentially: 

^'e * ^’eo ^ < 0 (3.7) 

where, is the random electron current density given in equation (3.3) 
and V is the probe potential. Potentials more negative than the plasma 
potential are electron retarding potentials; only those electrons with 
energies greater than the retarding potential can strike the electrode. 
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When ion currents can be ignored, it is a straightforward matter to 
determine the electron temperature by analyzing the electron component of 
the probe current as a function of probe potential for retarding potentials. 

If the active surface area of the probe is A then, in terms of the 
measured currents, equation (3.7) becomes 

- A 3 &xp{eV/kr ) (3.8) 

Taking the logarithm of both sides we have 

In i = ln(>l 3 ,) + eV/KT (3.9) 

e BO e 

Differentiating with respect to probe potential, V, gives 

^ any- JET C3.10) 

e 

whence 

-1 

(3.11) 

Therefore, if one makes a semilog plot of i versus V for electron retarding 
potentials the slope of the curve is inversely proportional to the electron 
temperature, from equation (3.9). This is the basis of electron-temperature 
measurement using retarding potential analysis. Figure 3.3. 

The problems associated with the method just described are experimental. 
The positive-ion current cannot be ignored. It is difficult to determine the 
precise passage of the probe through plasma potential, there being no obvious 
change in the electrode current. Since the method is valid only for retarding 
potentials, there is the possibility of applying the analysis to data taken 
beyond plasma potential and obtaining erroneous results (although perhaps 
not obviously in error). And finally, any type of system or natural noise 
must be minimized. 

To discuss experimental problems more fully and in preparation for 
Chapter 4, it is necessary to study the probe current -voltage relation more 
closely. This curve. Figure 3.4, is the probe characteristic. The electron 
current to a positively biassed electrode (or, the positive-ion current to 
negatively biassed electrode) depends on the probe geometry. For simple 
geometries exact expressions are available. For a planar electrode with 


T = 
e 


k d 

1 w ^ 
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Figure 3.3 (a) Theoretical electron current versus electrode 

voltage (spherical electrode); (b) Semi- logarithmic 
theoretical electron current versus electrode 
voltage (spherical electrode). 
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V(VOLTS) 


Figure 3.4 Total electrode current, including electron and 
positive-ion components. 
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dimensions which are large c<»npared to a Debye shielding length, the 
electron current for positive potentials is just the random current density 



(3.12) 


For a cylindrical electrode which is long relative to and of small cross 
section. 



(3.13) 


For a spherical electrode of small radius the currert for F > 0 is linear, as 
in Figure 3.3 



(3.14) 


The small spherical probe has an interesting characteristic: dj /dV is 

continuous through plasma potential. The plane and the small spherical 
electrodes are limiting geometries in that an electrode of any geometry will 
have a characteristic contained between these two characteristics. Equations 
(3. 12)-(3.14) are due to Mott-Snith and Langniuir [1926] and apply equally well 
to positive ions by changing the e subscript to +. The equations are then 
valid for positive ions for V < 0 (i.e. accelerating potentials). The 
principal assumptions here are that the plasma has a Maxwellian velocity 
distribution, that the plasma is at rest relative to the probe, and that no 
magnetic field is present. 

If we assume a spherical probe, we may readily calculate the magnitude 
of the retarding potential for which ion current becomes appreciable in 
comparison to electron current. Suppose the ion current is 10% of the 
electron current. Then 


A eV \ ^ . f eV \ 

y ' kT^ J ‘ ^eo J 


y < 0 


(3.15) 


The quantity |0.1 j /j | is about 17. 

that T a T - T, so that 
+ e 



For convenience, assume once more 


In 


(3.16) 
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The solution is 


et' 1 ft 

(3.17) 

or 


kT -4 

V • -1*8 ~ « -1.55 X 10^ T 

c 

(3.18) 


Clearly, at low temperatures the ion current becomes appreciable even for 
very small retarding potentials. For T • 300 K the ion current is 10% of 
the electron current for V • -47 mV. This is critical in that the 
uncertainty in locating the plasma potential is typically about SO mV. For 
T s 1000 K the ion current does not become 10% of the electron current until 
the probe is at -155 mV. 

Until about 1965 the following technique was used to handle the positive 
ion current. With the probe characteristic plotted as in Figure 3.4 the 
positive-ion component of the current (in the region between floating 
potential and plasma potential) was estimated by extrapolation of the curve 
obtained for large negative potentials, where it may be assumed the current 
is solely ion current. Subtracting the ion current from the total current 
yields the electron current in the retarding region which can then be 
analyzed by the standard sonilog plot, as described previously. In 
addition to making the data analysis quite tedious, two sources of error 
limit the accuracy of this technique to about ±100 K. The method involves 
subtracting two quantities of comparable magntidue to obtain a small 
residual. Errors in locating the plasma potential will be reflected as 
errors in temperature measurement. 

In 1967 a new technique for data reduction was developed [Smith et aX., 
1968] which resulted in accuracies of ±25 K under optimal conditions. The 
total current to a Langmuir probe of area A for electron-retarding potentials 
is i = where j ^ is the total current density, given by 

exp ^ + (3.19) 

e 

Again, the explicit form of as a function of V and is dependent upon 
the probe geometry. The derivative with respect to voltage is 


38 


•'« * „ •V' . 
nr ■ -0!- «P IT * ■3T 

^ B 

When dj^dV » ddJdV, equation (3.20) becomes 

■3T ■ "SJ" “P J 

Taking the logarithm of both sides of this equation, we have 

ln(ddJdy) - InCj- ^ e/kT) * eV/kT^ (3.22) 

which is similar to equation (3.9) 

A sonilog plot of dj^/dV will be a line with slope inversely 
proportional to electron temperature. As before, this is true only over a 
limited range in voltage from plasma potential down to the retarding 
potential where ion current becomes appreciable; this is precisely the 
constraint djJdV » djJdV. However, unlike before, the range of validity 
is easy to detect; a line may be fitted to the linear portion of the curve 
to immediately determine temperature. Figure 3.5 is a representative graph 
using this method of the data from Figure 3.4. 

The increase in accuracy using this technique comes from two sources. 
First, the range of validity of theory is operationally more clear. Second, 
Smith et al. [1968] differentiated the probe current on board the vehicle 
using a capacitor. The displacement current in the capacitor (which is 
proportional to dj /dV) was then measured using a logarithmic electrometer, 
itJ output being telemetered to ground on a conventional FM/FM telemetry 
link. Compressing the data in this manner prior to transmission makes the 
data more immune to noise and easier to process. The importance of this 
experimental arrangement must be stressed since it is ultimately noise which 
limits the accuracy of retarding potential analysis. It is this technique 
employing partial on-board data analysis that is used in the present research. 
The implementation of this experiment is described in the next chapter. 


(3.20) 


(3.21) 
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4. EXPERIMENTAL TE:HNIQUE 

4.1 Introduotion 

The Langmuir probe experiment together with several other, related 
experiments are integrated into the rocket payload. For the three flights 
that are the concern of this report the rocket is the spin>stablited Nike 
Apache, having a payload diameter of 6.S in (16. S cm). The configuration 
of this two-stage vehicle is shown in Figure 4.1. 

In addition to the experiments the payloads each include one or more 
batteries and their associated switching circuits and an FM/FM telemetry 
system having, typically, ten voltage-controlled subcarrier oscillators, 

(VCO's or SCO's), each representing one information channel. The 
frequencies of the various SCO's conform to IRIG standards. 

The channel assignments for Nike Apache 14.532, as an example, are 
listed in Table 4.1. The signals transmitted from this rocket are from the 
propagation experiment [Fillinger et al . , 1976] . the tip probe, including 
electron density fine structure [Klaue and Smith, 1978; ,'^ith and Klaus, 

1978] , and the boom probe and RF probe [Smith et al . , 1978] . The 
magnetometer, mounted transverse to the spin axis, gives aspect information 
and spin rate (used in the propagation experiment) . 

The telemetry channel assignments for Nike Apache 14.533 are given in 
Table 4.2. This payload contained two energetic-particle spectrometers, 
using solid-state detectors. One spectrometer (EES) was sensitive to 
electrons and protons; the other (EPS), by using a magnet in front of the 
detector, was sensitive only to protons. The payload of Nike Apache 14,534, 
launched earlier, was similar but included only one spectrometer, with no 
magnet . 

The emphasis of both flights was on airglow (391.4 nm) and energetic- 
particle flux. Stxne of the data relating to the energetic particles has been 
presented by Voaa and ^ith [1979]. The most significant result is the 
identification of the energetic particles as protons; they had earlier been 
assumed to be electrons. 

4.2 Probe Deeaription 

The instrumentation for the Langmuir probe experiment consists of three 
main parts: (1) the electrode and its supporting structure; (2) a programmed 

voltage applied to the electrode; and (3) a current -measuring device (electro- 
meter) which may also include logarithmic compression and differentiation. 



APACHE 


ADAPTER 
NIKE MOTOR 


Figure 4,1 Nike Apaches 14.532, 14.534, and 14.533. 



Table 4.1 Telemetry channel assignments for Nike Apache 14.532 


Sig”a 


RF probe, output 
Tip probe, fine structure 
R.-:eiver 41, modulation 
Receiver 42, modulation 
Tip probe, log output 
Boom probe, log output 
Boom probe, linear output 
RF probe, monitor 
Receiver 41, AGC 
Receiver 42, AGC 


Channel 

Number 

Center 

Frequency* 

(kHz) 

Information 

Bandwidth** 

(Ht) 

19 

93 

1395 

18 

70 

1050 

17 

52.5 

790 

16 

40 

600 

IS 

30 

450 

14 

22 

330 

13 

14.5 

220 

12 

10.5 

160 

11 

7.35 

110 

10 

5.4 

81 

9 

3.9 

59 


‘deviation ■ ±7.5% 
“modulation index ■ 5 
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Table 4.2 Telemetry channel assignments for Nike Apache 14.533. 


Channel 

Number 

Center 

Frequency^ 

(kHz) 

Information 

Bandwidth^^ 

(Hz) 

Signal 

21 

165 

2475 

Photometer, output 

20 

124 

1860 

EES, output #1 

19 

93 

1395 

EPS, output #1 

18 

70 

1050 

Tip probe, fine structure 

17 

52.5 

790 

Receiver #1, modulation 

16 

40 

600 

Receiver #2, modulation 

15 

30 

450 

Tip probe, log output 

14 

22 

330 

EES, output #2 

13 

14.5 

220 

EPS, output #2 

12 

10.5 

160 

EES, output #3 

11 

7.35 

110 

Photometer, filter temperature 

9 

3.9 

59 

Magnetometer 

7 

2.3 

35 

EPS, output #3 


Meviation = ±7.5 
••modulation index * 5 
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All three flights used the nose-tip electrode as shown in Figure 4.2. 

It has a base diameter of 1.08 in (2.75 cm) and a height of 1.50 in 
(3.81 cm). It is fabricated from a low-sulphur stainless steel (304 CRES) , 
intended for vacuum applications. To minimize outgassing during the 
flight it is baked in vacuum at a temperature between 900“ and 1000“ C 
for 30 minutes and then allowed to cool in the vacuum. 

The nose tip is the perferred location for a Langmuir probe because it 
minimizes, at least on ascent, the possibility of disturbance of the ambient 
plasma and of direct modulation of the probe current by the wake of the 
rocket as the rocket spins and precesses. On descent the rocket remains 
essentially upright to an altitude of about 100 km before beginning to tun* 
over. Some effects of rarefaction in the wake are noticeable in the descent 
data [Smith, 1969] . 

Electrons in the E region have only a very slight mobility across 
geomagnetic field lines. This is due to the rarity of collisions which 
would scatter electrons across field lines. Hence, electrons are constrained 
in their motion to follow helical trajectories along magnetic field lines 
with a radius of gyration (Larmor radius) of about 1 cm. At midlatitudes the 
rocket flight path places the rocket axis roughly parallel to the magnetic 
field. However, this is not always the case because of launch constraints 
and, in any event, the rocket will process during flight so that the cross- 
sectional area of the probe projected in the direction of the magnetic 
field will vary with time. To first order it is the projected area which 
is the effective probe area for electron sensing [Smith, 1969] . To 
minimize the sensitivity of the probe current to rocket attitude the tip 
electrode is ogival in shape. (An ogive is a solid formed by revolving 
an arc of a circle.) The height is chosen to yield a nearly constant 
cross-sectional area for forward angles. 

The tip probe is used for two types of data: electron temperature and 

electron concentration. To determine electron temperature the probe voltage 
is swept from 0 to 4.05 V with respect to the rocket body. This voltage 
range is sufficient to ensure the probe voltage has passed through plasma 
potential. Between probe sweeps is a period when the probe voltage is held 
at 4.05 V. In this constant-voltage mode the probe current is proportional 
to the ambient electron concentration. In conjunction with a radio- 
propagation experiment in the same payload (which serves to determine 


I 

I 


i 

i 

i 




} 



46 


the proportionality constant between probe current and electron concentra- 
tion) this mode of probe operation permits absolute measurement of the 
electron concentration with a spatial resolution of 1 m. Figure 4.3 is a 
graph of probe voltage versus time showing the periods over which temperature 
and concentration data are recorded. The tip-probe cycle period is 2 seconds 
with the swept-mode period being, here, 0.3 s (earlier 0.5 s was used). 

Since the rocket velocity is about 1.5 km s"^ at 100 km, there is one 
temperature measurement every 3 km. As the rocket travels upward the 
velocity continuously decreases resulting in an increased altitude resolution 
in the temperature measur^ents. 

In addition to a tip probe Nike Apache 14.532 (but not l*+.533 and 
14.534) also carried a boom-mounted probe as shown in Figure 4.4. 

The three-section boom is folded inside the payload at launch. A timer 
within the payload releases the boom cover plate when the rocket is at an 
altitude of about 60 km. The boom extends to a length of 80 cm by 
centrifugal force. At the end of the boom a spherical aluminum electrode, 

1 cm in diameter, and coated with Aquadag (colloidal graphite) is mounted 
on an insulator. The boom-mounted probe should have no magnetic aspect 
dependence since the electrode is spherical. The boom-probe on this flight 
was used only to measure electron tanperature. The boom-probe sweep 
program is also shown in Figure 4.3. The boom-probe voltage is swept from 
-1.35 to 4.05 V over a period of 0.5 s, with a repetition rate of 1.7 Hz. 

All current to the Langmuir probe must return to the plasma by way of 

the rocket body. It is desired that the ratio of rocket surface area to 

probe electrode area be large enough that even when the probe is positive 

relative to the plasma (and considerable electron current is being drawn 

by the probe) the floating potential of the rocket will not be significantly 

disturbed. In the E region the ratio of electron to positive-ion random 

current density is about 200. In order to restrict the voltage shift 

of the rocket relative to the plasma, it is usual to design rocket 

experiments such that the ratio of rocket surface area to probe area is at 

least an order of magnitude greater than J /j , or about 2000. The total 

eo *0 42 

area of an Apache second stage plus payload is about 2.6 x 10 cm while 

2 

the projected area of the ogive tip probe is 5.9 cm . This is an area ratio 
of 4400. To maintain this high area ratio and ensure a low impedance current 
path to the plasma, the Apache and instrumented payload are left unpainted. 


TIP-PROBE SWEEP PROGRAM 
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The full geometric area of the Apache and payload is used in this calculation 
(rather than the end-on area) because all the area is accessible to positive- 
ion current. This is due to the large Larmor radii of the positive ions 
('40 m for N* and O*). 

4.3 Probe Circuit and Operation 

The Langmuir probe circuit is shown in block form in Figure 4.5. 

With both relays in the positions shown a constant 4.05 V is applied to the 
probe. In this state the logarithmic electrometer monitors the probe current 
at constant potential. The output voltage is adjusted to fall in the range 
0 V to 5 V in order that it may be conveniently telemetered using conventional 
FM/FM telemetry equipment. This signal is used to obtain the electron- 
concentration profile as discussed earlier (Section 4.2). The output from 
the logarithmic electrometer also is coupled through a capacitor to an 
amplifier. As the rocket passes through small-scale plasma irregularities, 
or if a propagating plasma wave is present in the rocket environment, its 
spectrum will be contained in the output of the amplifier. 

After a period of 1.7 s the double-pole relay is energized for a period 
of 0.3 s by a trigger circuit which simultaneously enables the ramp generator. 
The probe voltage increases linearly with time over the 0.3 s from 0 to 
4.05 V. The probe current is monitored by the linear electrometer and applied 
to the logarithmic electrometer by way of a differentiating capacitor. In 
this mode the output from the logarithmic electrometer is \og{,di/dV) , the 
quantity which contains the electron temperature information. 

The above paragraphs refer to the tip-mounted probe. The boom-mounted 
probe in Nike Apache 14.532 does not use the constant-voltage mode. The 
ramp generator was enabled every 0.6 second and the double-pole relay was 
omitted. The boom probe and tip probe were not synchronized. 

All instruments are calibrated in flight by substituting a calibration 
resistor for the probe electrode. This is done by a calibration relay which 
is energized for two periods of about 5 seconds each. This relay is 
controlled by two barometric switches, the first being closed for altitudes 
greater than 40,000 feet (12 km) while the second is closed for altitudes 
less than 70,000 feet (21 km). With these switches in series the calibration 
relay is energized between 40,000 and 70,000 feet on both ascent and descent. 
The calibration resistor is typically 10^ Q. for daytime flights (14.532 and 
10^ n for nighttime flights (14.533, 14.534). In addition to the in-flight 


FLOATING GATT 16.2V 



Figure 4,5 Langmuir probe instrumentation block diagram 
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calibration an extensive calibration was performed prior to each launch. 

This consisted of sequentially coimecting resistors with values between 
10 4 

10 and 10 n between the electrode and ground. 

The electronic circuits of the probes, whether using the nose-tip or 
boom-mounted electrodes, are essentially identical. As noted above they 
differ principally in the waveform of the voltage applied to the electrode; 
this is determined by the timer circuit. 

The probe timer circuit is shown in Figure 4.6. The unijunction 
transistor Q1 and associated circuit sets the repetition rate of the sweep- 
fixed voltage sequence. Transistor Q2, which is normally conducting, 
energizing the relay in the sweep circuit, is cutoff for 20 ms after the 
period set by Ql. In setting up the times R2 is adjusted to give the 20 ms 
cutoff; then R8 is adjusted for 2 s (for the nose-tip probe) or 0.6 s (for 
the boom-mounted probe). These times are not critical and no attempt is 
made to temperature-compensate the circuit. 

The part of the circuit including transistors Q3, Q4, and Q5 is used 
to energize the mode relay for the duration of the sweep (0.3 or 0.5 s). 

The circuit including transistor Q6 is used to inhibit the sweep. Q6 
is in parallel with Q2 and, when conducting, prevents the sweep relay being 
de-energized and thus the sweep is not triggered. The sweep-inhibit circuit 
is controlled by a timer elsewhere in the payload. It is used in the nose-tip 
probe circuit to obtain a continuous (fixed-voltage) profile to an altitude 
of 120 km, after which the probe reverts to the alternating sweep/fixed- 
voltage mode. The sweep inhibit is turned on by the 70,000 ft (21 km) 
baroswitch so that the pre-flight and calibrated periods are also obtained 
in the fixed -voltage mode. 

The waveform generator which applies the voltage to the probe is shown 
in Figure 4.7. Since this circuit is inserted between the probe and the 
electrometer (current-to-voltage converter) it cannot be connected directly 
to ground. Accordingly it is powered by mercury batteries (12 cells of 1.35 V 
each). Until power is applied to the payload the sweep relay K1 is not 
energized, disconnecting the batteries from the circuit. When powev is applied 
to the circuit by energizing Kl, the output voltage (developed across R7 by 
the current is Q2) begins to rise at a linear rate determined by R1 + R3 and 
Cl. After 0.3 or 0.5 s, the value chosen for the sweep duration by adjusting 
R5, the voltage remains constant. It is adjusted by R6 to be 4.05 V 
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(chosen because it can also be obtained by three mercury batteries). 

When KI is de>energized (for 20 ns) by the timer circuit, capacitor Cl 
is discharged through resistor R4 (which limits the current). 

The relay K2. on being energized, disconnects the waveform generator 
from the electrode and substitutes a calibration resistor R8. As noted 
earlier this calibration period is controlled by baroswitches. 

The two arrangments of the floating-battery circuit are shown in 
Figure 4.8. The nose-tip probe usually has a 0 to 4.05 V sweep and the 
boom-mounted probe has a -1.35 to 4.05 V sweep. 

The PROBE IN terminal on the battery pack is connected to the input 
of the log-linear electrometer circuit, shown in Figure 4.9. Relays Kl and 
K2, when energized by the timer circuit, switch the input to the linear 
electrometer and connect the output of the linear electrometer through the 
differentiating capacitor C2 to the input of the log electrometer. When 
the sweep ends and the relays are de-energized the input is connected directly 
to the log electrometer, and the linear electrometer is disconnected from 
the circuit by opening the input. 

The linear electrometer uses a Fairchild F-418A operational amplifier 
with an RC feedback network. The resistor, R6, establishes the sensitivity: 
typical values are 1 Mn for daytime flights and 10 MO for nighttime flights. 
The parallel capacitor C3, has a value of 330 pF for daytime flights and 
33 pF for nighttime flights giving an upper 3 dB frequency for the 
electrometer circuit of about 500 Hz. The electrometer is followed by a 
741 op amp connected as a unity-gain inverting amplifier. 

The log electrometer uses a Keithley 302 operational amplifier with 
diodes in the feedback circuit to give the logarithmic characteristic. The 
voltage divider formed by R2 and R3 multiplies the voltage across the diodes 
by about 10 allowing a full-scale output greater than 5 V. 

The frequency response of the log electrometer varies with current. 

It is shown in Appendix I that the frequency response is better than 500 Hz 

-8 

for currents greatei than 2 x 10 A. 

When the probe is used strictly for electron-temperature measurement and 
the probe is swept repeatedly (without an extended period at fixed-voltage), 
as for the boom-mounted probe on Nike Apache 14.532, then the mode-switching 
relays are not needed and the circuits are correspondingly simplified. 





rometer. 













The main power supply for the payload is a re-chargeable battery 
consisting of 20 HR-1 cells. This 30 V supply (actually closer to 28 V) 
is reduced to 15 V (or other appropriate value) by voltage regulators in 
the individual experiments. A DC-to-DC converter is used to provide a 
-30 V supply for the payload. Again this is reduced by voltage regulators 
in the different circuits. 
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5. DATA REDUCTION 

5.1 Introduotion 

On beard the rocket the signal which carries the temperature information 
is developed as a voltage in the range 0 to 5 V. This voltage is applied to 
a voltage«controlled subcarrier oscillator (VCO) whose output frequency 
modulates a 1 W transmitter. The frequency of the transmitter is P band 
(about 240 MHz) or S band (about 2.2 GHz). On the ground this telemetry 
signal is received by a high-gain tracking antenna, demodulated, and 
recorded on an analog magnetic tape, together with timing and other signals. 
As the analog magnetic tape is being produced the demodulated signal is 
also applied to frequency-to-voltage converters (discriminators) whose 
outputs reproduce the signals applied to the VCOs in the rocket. These 
varying dc signals are digitized and recorded on another magnetic tape. 

At any time after a launch has occurred the analog tape may be played back, 
discriminated, and digitized to produce additional digital recordings. 

The digitizer operates at a rate of 5 kHz and is preceded by a five- 
channel commutator. Every fifth data sample to be digitized and recorded 
is electron temperature data; the other four data samples are taken from 
other information channels. This results in a digitization rate of 1 kHz 
per data channel. 

Before digital ccanputers were widely available, all rocket data 
reduction was accomplished by hand. The recovered voltage from the telemetry 
discriminator was recorded on a high-speed stripchart recorder. For each 
probe aweep a temperature could be determined by measuring the slope of the 
leg {dUdv) signal using a straight edge and protractor. Figure 5.1 is a 
stripchart record demonstrating this technique. Data reduction was 
straightforward, although tedious, and quite satisfactory under optimal 
C( iditions. 

The initial desire to computerize the data reduction process was to 
save time and reduce errors. After software was created to handle the 
digital magnetic tapes other advantages became apparent, foremost being the 
capability to incorporate simple noise-reduction algorithms. This permitted 
the recovery of temperature information from data recorded under less than 
optimal conditions. 

The software must be able to perform several tasks: 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 5.1 Stripchart record demonstrating manual technique of recovering 
electron temperature information. 
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1. The program must efficiently scan the digital magnetic tape 
and recognize probe sweep data; 

2. The probe voltage sweeps linearly from 0 to 4.05 V. Only 
data generated while the probe is within (roughly) 100 mV 
of plasma potential are meaningful for temperature 
determination. Other data must be ignored; 

3. Once the proper data have been located the temperature 
determination must be made. This entails fitting a line 
to a limited number of data points. Data which show 
evidence of systematic error must be corrected at this stage 
or eliminated; and 

4. With rocket trajectory information supplied by NASA the 
program must assign an altitude to each temperature 
measurement . 

After the temperature data have been generated they should be printed in 
tabular form and, ultimately, plotted as a graph of temperature versus altitude. 
5.2 Data Reduction Software 

The program to obtain electron temperature is ELTEMP. Before proceeding 
with a detailed description of the function of each subroutine in ELTEMP it 
is necessary to describe the organization of information on the digital tapes. 
Each digital tape consists of one file containing records of various lengths. 

At the beginning of each tape are four header records, each containing 45 
words of identifying information. Each word contains 16 bits, the equivalent 
of an IBM 360 half-length integer. Alpha-muneric information is encoded 
using the IBM Corporation EBCIDIC convention. 

Following these header records are five calibration records, each 
containing 1005 words. These records are the digital levels corresponding 
to the five analog telemetry levels, 0 V, 1.25 V, 2.50 V, 3.75 V, and 5 V. 

The remaining recoras are all data records, containing 2008 words each. 

This organization is demonstrated in Figure 5.2. 

The data record structure is displayed in Figure 5.3. Each of these 
records contains five channels of data ninning cyclically through the 2008 
words, beginning with word 6 and ending with word 2005. Words 1 through 5 
of each record are not used. Words 2006, 2007, 2008 contain binary-coded 
decimal time as explained in Table 5.1. The time recorded in these three 
words is the time of word 6 in the following record. Note that there is a 
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FILE 

1 

1 

2 

3 

4 

5 

e 

7 

8 

9 

10 

II 

12 

13 


^ ® ^ ^ 


END OF 
FILE 


^ Four header records of length 45 words containing tape identification, 
^ e.g., AMQl NIKE APACHE 14.532, 

Five calibration records of length 1005 words. These records are the 
(§) digital levels corresponding to the five analog telemetry levels 
LBE(-7.50%), LBH(-3,75%), BAND CENTER, UBH(+3.75%), and UBE(+7.50%). 

© Data records of length 2008 words. 

Figure 5,2 Organization of data on a digital tape. 



















Table 5.1 NASA format for BCD time encoding 


Word 

2006 bits 



1 - 4 

tenths of milliseconds 


5 - 8 

one of milliseconds 


9-12 

tens of milliseconds 


13 - 16 

hundreds of milliseconds 

Word 

2007 bits 



1-4 

ones of seconds 


5-7 

tens of seconds 


9-12 

ones of minutes 


13 - 15 

tens of minutes 

Word 

2008 bits 



1 - 4 

ones of hours 


5-6 

tens of hours 


7-10 

ones of days 


11 - 14 

tens of days 


15 - 16 

hiindreds of days 
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constant increment of time between any two adjacent words in the same channel 
of 1 ms, even between word 2005 of one record and word 10 in the following 
record. Since each data record contains 2000 words of data digitized at 5 kHz, 
one record corresponds to 0.4 second of data. This is 400 words of 
temperature data and 1600 words frcmi the other four channels. 

Program ELTEMP consists of the nine major FORTRAN subroutines listed 
in Figure S.4 allowing it to operate on the control Data Corporation CYBER 
computer. The main program is simply a calling program for the many 
subroutines. The entire program is listed in Appendix II. 

The tape information is read in subroutine LOAD and the logarithmic 
electrometer information, log (dt/iy), is transferred to the array L0G(Af,/l/). 

The data in the other channels are discarded. This array is dimensioned 
M - 3, /? = 400 permitting three data records of 400 points to be stored. 

Each time LOAD is called one data record is read from the tape and stored 
in LOG with M equal to 1, 2, or 3. LOG is a circular file in that it is 
filled sequentially from the magnetic tape as shown in Figure 5.5. 

Called from within LOAD are the subroutines CALTIM and BINARY. BINARY 
converts the words 2006, 2007, and 2008 into binary representation. CALTIM 
takes these three binary numbers and calculates the time using the timing 
code in Table 5.1. This is the time (UT) in days, hours, minutes, seconds, 
and fractional seconds. From this is subtracted the time of launch to give 
the time after launch (in seconds). This is the tape record time which is 
stored in RECTIM(M) upon returning to LOAD. Also, the tape record number 
is stored as RECNUM(M). 

The main program then calls subroutine TRSHLD which searches through 
the data is array LOG for the beginning of a probe sweep. The signal 
characteristic which indicates the beginning of a probe sweep is the sharp 
transition (<5 ms) from a high signal value KMAX to less than a small signal 
KMIN for 10 ms or longer. The values l^iAX and KMIN are input variables 
since they vary from flight to flight depending on the digitization constant 
(digital units/analog volt) used at the ground telemetry station. 

In practice it is a simple matter to have the contents of a portion of 
data tape printed, observe a transition at the beginning of a probe sweep, 
and choose appropriate values for KMAX and KMIN. A useful program for this 
purpose is NASA listed in Appendix III. Ttic signal must stay low for more 
than 10 ms (10 points) to remove the possibility of a single low bad point 















Figure 5.5 Circular data procersing file. Logarithmic 
electrometer data are shifted into array 
LOG 400 points at a time. Associated with 
each array LOG is RECNUM, the magnetic tape 
record numbci , and RliCl IM, the record time 
in seconds after launch. 
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triggering the program falsely. The testing in TRSHLD is sufficiently 
general to operate across an array boundary, say between LOG (2, 400) and 
L0G(3,1), if the transition should by chance occur there. 

If TRSHLD does not find a sharp transition the program control returns to 
the main program. Here the central index M is incremented by 1. If W » 3, it 
is incremented to M = 1. Subroutine LOAD enters another record of probe data 
in array LOG. CALTIM calculates another time and TRSHLD is called once more. 

When TRSHLD does locate the beginning of a probe sweep program control 
passes to CRUNCH, the most complex of the subroutines. The purpose of 
CRUNCH is to determine the slope of the log idijdv) data about plasma potential 
(in the retarding potential region), take the reciprocal, and multiply by a 
calibration constant to recover the temperature in degrees Kelvin. The first 
action CRUNCH takes is to advance some number of points from the transition 
to the center of the probe sweep. The number of points advanced is an input 
parameter called OFFSET. During the rocket flight the vehicle acquires 
floating potential a few volts belo” plasma potential. When the probe sweep 
begins (at zero volts with respect to the vehicle) there is very little 
electron current because the probe is well below plasma potential. There is 
no reason to examine this portion of the probe characteristic. The OFFSET 
parameter locates the region of interest closer to the point when the probe 
passes through plasma potential. 

The next 60 data points are that portion of the probe characteristic 
containing temperature information. As recorded on the digital tape this 
may bridge two data records and hence be transferred to two arrays by LOAD, 
say L0G(2,-) and L0G(3,-). For computational convenience CRUNCH transfers 
these points to one 60-element array LLOG. The calculating window was 
chosen to be sixty points in length; this much data permits one to get an 
intuitive feel for the probe characteristic. Also, using fewer data would 
cause the choice of the OFFSET parameter to be quite critical. As the rocket 
ascends it usually encounters higher temperatures causing the floating 
potential to be lower. This means the probe passes through plasma potential 
further into the probe sweep. For Nike Apache 14.532 and 14.534 the 60- 
point data window was wide enough that the OFFSET parameter did not need to 
be adjusted during the flight. 

With the raw data now in array LLOG subroutine CRUNCH proceeds with a 
five-point line fit, first with points 1 through 5, then 2 through 6, and 
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SO on. nils genorates S6 slopas tdiich ars stored in array SLOPEZ. The line 
fit is a least-squares fit and incorporates the following noise-reduction 
technique. Five points are used in the line fit. Then each point is checked 
to see how far it lies from the best-fit line. If this distance is less than 
the number SPREAD (usually 2S0 units), CRUNOI will proceed with the next five- 
point line fit. If the distance is greater than SPREAD, that point is thrown 
out and the line fit done again with four points. Again the points are 
checked to see how far off the line they lie. If any of the four points fall 
beyond SPREAD from the best*fit line, the slope for those points is set to 1 
and CRUNCH proceeds with the next five poi.nt line fit. 

When the line fitting is complete for a given probe sweep there are 
56 slopes in the array SLOPEZ. A running average of the slopes (5 points 
at a time) in array SLOPEZ produces 52 smoothed slopes which are stored in 
SLOPEZ again. The reciprocal of each slope is multiplied by a calibration 
constant to yield the numbers in array TE, one of which is the electron 
temperature in degrees Kelvin. Subroutine MINI makes this choice by 
choosing the lowest temperature as representative. Although simple, this 
choice is completely satisfactory. Figure 5.6 shows the contents of arrays 
LLOG(in digital units) and TE (in degrees Kelvin) for a probe sweep during 
the flight of Nike Apache 14.532. As can be seen there is an extended region 
of 8 ms where the slope is constant and the temperature is well defined. 

Each probe sweep is printed in this fashion to verify the quality of the data 
and as a diagnostic when one temperature in a profile is far from the mean. 

The quality of the data is actually better than Figure 5.6 appears for the 
figure was produced on a cnnputer drum printer resulting in a quantitzation of 
the data. The printer plot is produced by subroutine PLOT. The time of the 
temperature measurement is determined by subroutine ITIME by interpolating 
between the times given at the end of each tape record. 

Each call to CRUNCH results in a time and temperature being returned 
to the main program. This continues until the tape recorder number matches 
the stop number given by the user. Subroutine TRAJ performs a parabolic 
interpolation using the NASA trajectory data to assign an altitude to each 
time/temperature pair. Finally a table of time, altitude, and temperature 
is printed. 



ttc(i aRainst word mimbcr. The iriniiPur! value in the center of the 
ure of the Tl data is the electron tenperature (K ) during the prohe 
cp cycle. Note that one word corresponds to a one millisecond 
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S.3 Calibration 

In the retarding potential region for an electrode of any general shape 
the electron current density to the probe surface is 

0 

where J is the random electron current density to the probe surface 
00 

with the probe at plasma potential, k is Boltzmann's constant, e is the 
electronic charge, and T is the electron temperature. The voltage is 
negative voltage measured, strictly speaking, from the plasma potential. 

In the experiment the probe voltage is measured from the Apache rocket body 
which floats about 1 V negative with respect to plasma potential. As will 
be seen this offset is unimportant so long as it does not vary during a 
probe sweep. 

The probe has an effective area for current collection. A, thus the 
total electron current to the probe is 

t exp(j^) ^5.2) 

0 


The voltage applied to the probe is generated by a ramp generator 
(Figure 4.7). Let the voltage sweep rate be and the time from the 
beginning of the sweep be t (the voltage is assumed to be 0 V at t«0; 
the actual value is inconsequential) then 

^ J’eo ®*P 

This current is converted to a voltage by a linear electrometer with an 
effective (feedback) resistance /? so that the output voltage of the linear 
electrometer is 


ev^t 


v(t) m RA exp 


(k'/' 




(5.4) 


From the linear electrometer the signal passes through a capacitor which 
time-differentiates the signal. For a capacitance C the displacement 
current is 


i(t) - cm 

The current through the capacitor is converted to a voltage by a second 
electrometer, this one employing a diitde as a feedback element to provide ;i 





logarithmic scale. For a silicon diode with a forward current' of less than 
5 lO"^ A the diode current is an exponential function of voltage. The 
transfer function of the log electrometer is 




(5.6) 


.-11 


in which has a value about 5 x 10 A, is the output voltage of the 
log electrometer and u is the e-folding voltage, typically 0.4 V. 


Thus 


u = u In 
o e 




(5.7) 


The input current is the capacitor displacement current so that the output 
voltage of the log electrometer is 


y (t) = y In 
c e 



(5.8) 


or 


“o'*’ ' ■■’e 


ft)l 


ey t 
8 

kT 


(5.9) 


This voltage is telemetered to ground. 

At the ground telemetry station a voltage is recovered from the 
discriminator. It is related to the voltage on-board the rocket by 
V = Gv^ + H, where is a scaling factor and ^ is a constant offset voltage. 
This voltage is carried over a coaxial cable to a digitizer operating at 
1 kHz per data channel. The digitizer has a digitization constant p (digital 
units/analog volt) and a positive offset r. The number recorded on the 
digital tape is thus 


u = p(Gv^ * H) * r> (5.10) 

where y^ is the voltage given in (5.9). 

To recover the electron temperature information from this signal the 
fixSt step taken is tc calculate the digital slope of the signal, i.e., the 
number of digital units increase per digital word. This eliminates the 
digitization offset r and the telemetry offset voltage H. 

Au = vGAv = pG 
o ' 

where q is the interval between digital samples (i.e. 1 ms). 



(5.11J 
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As explained in Section 5.1, tu is obtained in subroutine CRUNCH by a 
S'point line fit rather than by differencing adjacent data points. Dividing 
by p, the telemetry constant, yields analog volts per digital word. Dividing 
this by q, the digitization time increment, gives dojdt (volts per second) 


Am 

P<? 



(5.12) 


From equation (5.9) 



(5.13) 


In this equation C, F, I and are constants. The effective area of 
the probe. A, may change slightly during flight due to changes in the probe 
surface. However, the probe sweep rate is about 70 ms V"^ and only 200 mV 
of data are required in the retarding potential region. This is about 14 ms 
of probe data for which the assumption that A is constant is justified. The 
random current density j also changes during the flight, since it is a 
function of electron concentration. But again, with the vehicle traveling 
at less than 1.4 km s above 90 km there would have to be a large gradient 

in n in a distance less than 20 meters to significantly change d. The 

6 , &0 

only ionospheric structure capable of generating such a gradient is the 
sporadic-E layer. In all other cases can be regarded as a constant for 
the duration of the probe sweep. The only other variable is the electron 
temperature. There has recently been evidence indicating sharp temperature 
gradients within sporadic-E layers [Szuazczewvoz and. Holmes ^ 1977] and we 
must rule out the use of this technique in such circumstances. However, 
it is safe to say is constant over the duration of a probe sweep in all 
other situations. The time derivative of the log term is thus zero for 
all of the E region except, perhaps, in the sporadic-E layers. 

Combining equations (5.12) and (5.13) we have 


Am 

P<7 


Gv ev 
e 8 

kT 


(5.14) 


whence 


T = 
e 


kAu 


(5.15) 
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The notation used is as follows: e is the electronic charge (1.602 x 10 C), 

y is the probe voltage sweep rate (V s ^), p is the digitization constant 
® -1 

(digital units V ), q is the digitization time increment (s), k is 

-23 -1 

Boltzmann's constant (1 . 380 x lo J K ) and u is the measured 

e 

logarithmic e-folding voltage. The telemetry constant G is, in practice, 
combined with p . 

In subroutine CRUNCH the digital slope is calculated, its reciprocal 
taken and multiplied by CALIB, the calibration constant. From (5.15) it is 
seen CALIB is numerically equal to Gv^eo^q/k. The probe voltage sweep 
rate and the elecr^ometer e-folding voltage are measured while the 

payload is on the ground. The digital time increment q is 1 ms. Records 5 
through 9 on each digital magnetic tape contain the information to determine 
Gp. (The two quantities need not be measured separately.) These records are 
the telemetry calibration records: record 5 is the digital level for 0 V, 

record 6 is the digital level for 1.25 V, and so on, with record 9 being the 
digital level for 5 V. To determine Gp for channel 1 one can average all 
200 of the channel 1 values in record 9. Also average all the channel 1 
values in record 5. The difference divided by 5 V is Gp for channel 1. 
Program CALIB, listed in Appendix IV, was designed for this purpose. It 
prints the average channel values for records 5, 6, 7, 8, and 9. Examples 
of the output are shown. 

The e-folding voltage is determined from calibrations performed prior 

to launch (see Section 4.3). The probe is connected to the vehicle ground 

4 10 

by way of a series of test resistors ranging from 10 0 to 10 0. For each 
resistor, place in the circuit the resultant output voltage from the 

logarithmic electrometer is, 

CEV 

y = In -5 1 (5.16) 

which may easily be seen by going back to equation (S.2) and starting with 

i - VIE A semi-log plot of V against E , ^ has a slope v allowing 

oaL (iCii' c ' 

y to be obtained. 
e 

During the in-flight calibration periods a calibration resistor is 
substituted for the probe electrode. This gives a one-point check of the 
calibration. Experience in many flights has shown that this is sufficient 
verification of the system. 


74 


6. ERROR ANALYSIS 

This chapter discusses noise encountered in rocket experiments with 
Langmuir probes and the specific efforts taken to reduce the effects of 
noise in the data analysis. Five types of noise (error) have been 
identified and are considered in the following sections. 

1. Statistical scatter among data points. 

2. Power line noise. 

3. Systematic error in data handling and ref ixtion (low-pass 
filtering) . 

4. Probe surface contaiiiination (work function or contact 
potential) . 

6.1 Statistical Soattev 

The simplest type of noise encountered is statistical scatter on the 
data points. The source of this noise is probably within the payloads' 
instrumentation, including the experiment itself and the telemetry system. 

The scatter is of such a magnitude that it is not satisfactory to take the 
difference of adjacent points to evaluate the slope of log {dijdo). 

Therefore subroutine CRUNCH determines the slope by a five-point least-squares 
fit described earlier. This technique appears to be completely sati.sfactory 
under normal conditions. 

Severe scatter was observed on five rocket flights during 1974 and 1975 
in which there was a partial failure of the P-band telemetry transmitters 
shortly after launch. During the first few seconds of flight it is usual 
for rockets to acquire a static charge due to rocket motor firing and, to 
a lesser extent, triboelectricity. This charge causes voltage on the 
order of 2 kV to develop on antennas or probes extending from the vehicle. 

The rocket flights on which transmitter failures occurred used telemetry 
transmitters without a dc shunt in the output stage. The high voltages 
encountered during the first five seconds of flight were sufficient to 
damage the transmitter output stage so that the rf output was reduced to 
a few milliwatts. This has been experimentally verified by Gilchrist 
[in Edwar>ds, 1977] in a specially instrumented payload. 

Nike Apache 14.532 was affected by such a failure. The telcnctr) signal 
was poor throughout the flight making the recovered signal very noisy. The 
boom probe log electrometer data was usable until 173 seconds after launch 
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(182 km) and then became much worse, as illustrated in Fi^'ure 6,1. Subroutine 
CRUNCH, in doing the line fitting to determine slope, th’-ew out many data 
points as lying too far from the best -fit line. Even so, the data from each 
probe sweep had to be viewed by hand (using the plots generated by subroutine 
PLOT) to be assured the recovered electron temperatures were reasonable. 

The telemetry transmitter has been redesigned and flown successfully on 
following flights. Also, it is worth noting that the charge acquired by the 
vechile early in flight is dissipated well below 50 km and in no way affects 
the Langmuir probe. 

6.2 Power Line Noise 

Vfhile using program EI.TEMP on magnetic tape data from Nike Apache 14.534 
it became clear the data were contaminated with a periodic signal not related 
to rocket spin or precession. A 16 or 17 ms variation can be seen in Figure 
6.2 in the plot of reciprocal slope, although not in the data because of 
the poor resolution of the printer plot. In the telemetry building at 
Wallops Island Station the data lines are coaxial cables running from the 
telemetry discriminators through a patch panel to an A-to-D converter in the 
computer room, a distance of about 10 meters. It appears that somewhere 
along this cable run there was sufficient capacitive coupling to 6 O-H 2 power 
lines (t = 16.7 ms) to cause modulation of the analog data signals. To 
establish this, the telemetry step calibration at the beginning of two data 
tapes were plotted on an expanded scale. In Figure 6.3 are the telemetry 
step calibrations for Nike Apache 14.534 for the log electrometer channel, 
low pass filtered at 81 Hz. Figure 6.4 shows the telemetry step calibrations 
for the same channel on another tape, low pass filtered at 450 Hz. The 
traces are the digital values representing 0 V, 1.25 V, 2.50 V, 3.75 V, and 
5.00 V and should be perfectly constant. Instead, a strong CO-Hz component 
is seen on both tapes. In Figure 6.3 the peak-to-peak deviation is about 
24 digital units (32 mV). In Figure 6.4 the 60-Hz component is present at 
about the same mangitude, as well as a noticeable 120-Hz harmonic. 

Figure 6.5 is a plot of a portion of the in-flight (ascent) cal ibration 
for Nike Apache 14.554 at 23 seconds after launch. In the lower trace the 
maximum peak-to-pcak displacement is 52 digital units (67.5 mV). The 
seriousness of this can be understood by examining figure 6.2 again. The 
plc.r of reciprocal slope may be read as temperature in degrees Kelvin on the 
right-hand scale. By chance, just where the minimum temperature should have 



Reproduction of plot from record 479 from Nike Apache 14.532 showing noise 
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TELEMETRY' UNITS 




TELEMETRY UNITS 





Figure 6.5 A portion of the inflight (ascent) calibration for Nike Apache 14.534. The trace 
should, ideally, be completely constant. The modulation of the signal, 67.5 mV 
p-p, occurred at the telemetry station through capacitive coupling between the 
data cable and 60-Hz power lines. The signal is shown using two low-pass filters 
(450 Hz and 81 Hz) in the telemetry play-back system. 
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occurred there is instead a bump, attributable to contamination of the signal 
at 60 Hz. (Compare with Figure 6.5, an example of a probe sweep without 
60 Hz contamination.) It is clear, depending on the phase of the 60*Hz 
signal during the probe sweep, the recovered temperature can be elevated or 
depressed by as much as several hundred degree. 

Tests conducted at the Wallops Island telemetry facility showed that, 
with the discriminator in the telemetry room and the digitizer in the 
cranputer room there was a 30-mV p-p 60-Hz signal on the cable at the 
discriminator input. When the discriminator and low-pass filter were 
relocated adjacent to the digitizer and a very '.ior' interconnecting cable 
was used, it was found that this reduced the hum pickup to 10-mV p-p. In 
both cases the telemetry discriminator was receiving a band-center signal 
and was followed by a 59-Hz low-pass filter. Considering the short length 
of cab’* in the second case, it is likely that the residual hum 
originated within the discriminator or digitizer. 

Rather than redigitize the tapes at the 10-mV noise level, the decision 
was made to use the given tapes with software designed to tolerate (or 
remove) 60-Hz hum from the data, a process generally referred to is 
"humbucking". This approach has the advantage that it might be used on 
tapes from other launches which have been recorded with hum (such as Nike 
Tomahawk 18.1006 and 18.1008). Two humbucking algorithms were developed 
and tested on the Nike Apache 14.534 tapes. 

The first technique used the idea that when the probe is strongly 
negative it should draw only positive ion-current. Furthermore, since the 
ion-rs.ndom-cuiTent density is about 2U0 times smaller than the electron- 
random-current density, this portion of the characteristic may be taken as a 
rough baseline ^>r electron current and as a 60-Hz hum reference. The 
following algorithm was used: subtract from each data point the point which 

occurred 50 ms earlier. This should yield a curve without hum since 50 ms 
is the duration of precisely 3 cycles at 60 Hz. And although the curve would 
be offset by a constant (the positive-ion current) the slope would be pre- 
served. The improvement in data quality is shown in Figures 6.6 (a) and (b) . 

The above technique was used with confidence until an attempt was made 
later to quantify the positive-ion current to the probe for potentials less 
than plasma potential (i.e., at the beginning of the proiie sweep) It was 
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Figure 6.6 A comparison of two techniques for removing 60 !lz hum 
from probe sweep data. (a) original data with a 60-lU 
component visible in the temperature curve (b) 

and (c) same data after processing with 50 ms and 200 
ms, respectively, humbucking techniques descrii)ed in 
the text; (c) appears to be completely free of contam- 
ination. 


ELECTRON TEMPERATURE (K) -ELECTRON TEMPERATURE IK) —ELECTRON TEMPERATURE CK) 
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then discovered that the positive>ion current had a slight negative slope 
until the electron current became dominant further into the sweep. This 
slight negative slope, if subtracted from the electron*retarding portion 
of the characteristic, could artifically increase the slope there. The 
result would be to lower the apparent temperature. 

In light of this a second humbucking technique was developed. Before 
the beginning of each probe sweep is a period when the probe is held at 
a positive potential (♦4.05 V). This signal is constant over periods of 
time of interest here (>60 ms) except for 60 Hz hum. Therefore, this 
signal may be subtracted from the electron-retarding region to yield data 
relatively free of 60-Hz hum while preserving the slope. This was done 
in Figure 6.6 (c) where a 200-ms delay was used. The improvement in 
data quality is striking. The 200-ms humbucking routine resulted in 
temperatures averaging 135 K higher than the 50-ms routine. Note that 
in Figure 6.6 (c) is 125 K greater than that in Figure 6.6 (b) . 

The 200-ms technique worked well and, it is felt, accurately recovered 
the electron-temperature information. The computer code for this algorithm 
is documented in Appendix II. 

A simpler technique was used with the data from Nike Apache 14.533. 

There the electron temperature was so high that the constant slope persisted 
for over 80 ms. Thus it was reasonable to average 17 points at a time to 
produce data free of 60-Hz hum. The coding for this was incorporated 
within subroutine CRUNCH and is documented in Appendix V, 

6.3 LoD-Paes Filtering 

The signal from the Langmuir ,>robe is low-pass filtered scvcr.al times 
prior to and during data reduction. Low-pass filtering at intermediate 
stag's restricts high-frequency noise (>4S0 Hz) and permits the time 
der'vativcs necessary for temperature determination to be made. If the 
frequency response of any filtering stage is so low as to reduce the slope 
of the log(cft/<iu) signal, the result will be an apparent increase in the 
electron temperature. This effect might be important in the lower portion 
of the F region {--lOO km) where the electron temperature is about ."On k. 
Figure 6.7 indicates the various stages in which iow-pass filtering occurs. 

It is important to understand the distortion of the probe signal that 
occurs in the electrometer circuitry and telemetry systoiii. For the followirii’ 
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development refer to Figure 6.8, a simplified schematic of the probe 
electrometer circuit. For this discussion we assume that the current from 
the probe is strictly electron current, given by 



( 6 . 2 ) 


where A = effective probe area 

j » electron random current density 
eo 

Vjf » rocket floating potential 

V = probe voltage sweep rate 

e - electronic charge 

k - Boltzmann's constant 

» electron temperature 

Let a * (ev /kT ). The probe current flows into the first linear electrometer, 
which has feedback elements Rj and C^. The resultant output voltage is: 




A 3^^ exp 


e(v t - 
' 8 

lA. 

kT 


e 



Cj(o + 1/RjCp 




( 6 . 3 ) 


The electrometer is actually an inverting op amp, although the sign is not 
displayed explicitly. This voltage causes a current ij to pass through 
Rq into the second linear electrometer. The output voltage from it is: 







*^0^2 ^R^ " 
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typical values for the constants in the above expression are 


10 KO 


5.1 MO 


R2 = 10 KQ 


156717 K^S"^ 


68 pP 


C 2 » .01 mP 


1/RjCj 

I/R2C2 


2884 s 


10000 s 


Substitution of these values in equation (6.3) shows that the second and 

third terms in the RHS are negligible compared with the first term. 

The probe is 0 V for 20 ms before the sweep begins. This permits 

transients to die away and establishes initial conditions Vj(0") and 1 ^ 2 ( 0 "). 

During this period the only current to the probe is a very small constant 

-9 _ 

positive-ion current of about 5 x 10 A. This sets Vj(0 ) and 5 

25 mV. A J is about 1 lO"^ A. If T - 400 K, y„ = 13.5 V s"\ and 
€0 e 8 

Uy. 3 l.S V, say, equation (6.3) becomes 

V 2 (t:‘^) - -0.43 exp[29(13.St - l.S)] + 0.035 exp(-2884t) 

♦ 0.025 exp(-10"^t) (6.4) 


The first term is the response to the electron current and becomes 

-8 

significant at about t ■ 70 ms (V 2 * -4 x 10 V). The other terms are less 

than this at t > 5 ms and are completely insignificant during the time 
period of interest here. 

Voltage V 2 is differentiated by resistor Rj in series with capacitor 
Cj. The purpose of the differentiation is to ensure that the current ^2 
the logarithmic electrometer does not pass through zero during the sweep. 
(This would otherwise occur for, as mentioned earlier, at the beginning of 
the probe sweep there is some positive ion current.) Now, ignoring the 
transient terms in equation (6.3), current t 2 is 


A 0 a exp 
eo 


- v) 


i2(*) 


‘^o‘^3^1^“ RjCj ^ ^2^“ * R2C2 ^ ^ R3C3 

Vj(O') 


(6.5) 


where V_(0 ) is the voltage across C_ at t = 0. Again, the transient term 

^ -1 -1 
is unimportant because ® [330 Kfi x o.OOl pF] = 3030 s 
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Current ^ enters the logarithmic electrometer formed by an operational 
amplifier with diode feedback elements. The output voltage is a ramp 
containing the temperature information, with another texm: 


Vj(t) 


«(v - ■’/) 

y 






( 6 . 6 ) 


* V In 
e 


* R2C2^ RjCj^ 


where is the electrometer e-folding voltage and Iq is the reverse 
saturation current of the diode used as a feedback element. The second term 
is temperature dependent (recall that o = ev /kT ) and will be considered 
later. The first term is a ramp with slope inversely proportional to T^. 

Vj is telemetered to ground, low-pass filtered, digitized, and recorded 
on digital magnetic tape for later analysis. 

The error introduced by the low-pass filter may be significant, depending 
on T and u (the filter cutoff frequency). The filter phase and amplitude 
response is given in Figure 6.9 in a plot provided by the manufacturer, Data- 
Control Systems, Inc.* DCS refers to the filter as a third-order "linear phase" 
output filter, a term used to signify a Bessel response filter. The transfer 
function for a third-order Bessel filter, normalized for a dc delay of Is, is** 

(6.7) 

H(s) 


15 


15 


n=3 


s^ + 6s^ + 15s + 15 


(s + 3.67782s + 6.45944) (s * 2.32219) 


This may be renormalized for a cutoff frequency wal by substituting 1.753s in 
place of s in equation (6.7). This yields the actual transfer function for 
the response in Figure 6.9. 

2.78449 


H(s) 

«*1 


( 6 . 8 ) 


(s + 2.09801s + 2. 10199) (s + 1.32469) 


The response Y(s) from such a filter (with cutoff frequency w) upon the 
application of an input signal at is, 


*Data-Control Systems, Inc., P.O. Box 860, Commerce Drive, Danbury 
Connecticut. 

**see, for example. Principles of Active Network Synethesis and Design , 
Gobind Daryanani, John Wiley and Sons, 1976, p.g. 119-120. 



Figure 6.9 


Normalized frequency response characteristics of linear phase 
I'Z'lTotu" '^-^-Control Systems telemetry discrLiH^tC 


PHASE SHIFT IN DEGREES 
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Y(S) 


2.78449o)‘ 


(s^ + 2.09801US + 2. 1019901^) (s ♦ 1.32469o)) 


(6.9) 


This may be written as a partial fraction expansion (6.10) 

_ -1.753 a/iii . a . f. 28046 a/hi )s ♦ 1.53907a . 1.47254 a/ui 

« _________ 4. + ♦ ' 

s s s ♦ 2.09801US ♦ 2.10199a) s + 1.32469u) 

The time domain response can then be determined from a table of Laplace 
transforms. For an input signal at u(t) the output is (omitting «(*)), 

Y(t) » -1.753 o/o) + ot + 1.275511 a/u (6.11) 

COS(-l. 34902 + 1.00079o)t) + 1.47254 a/o) g“l-32469a)t 
The input ramp comes through the low-pass filter with slope unchanged, but 
with three additional terms: a dc offset, an overdamped oscillation, and a 

decaying exponential. The magnitude of these terms is proportional to a/oi. 

It is the slope of the response (albiet modified) which contains the 
temperature information. Taking the time derivative of equation (6.11) yields. 

-1.04900o)t 


Y(t) » a[l - 1.33760 e 
. 1.27612 
- 1.95066 


cos(-1..32t902 + 1.00079o)t) 
sin (-1.34902 + 1.00079o)t) 

( 6 . 12 ) 

Here a is the desired slope, the remaining terms being error terms. 

It is interesting that the error terms depend only on the product oit and 
display no dependence on temperature. Measuring the signal slope at t = 15 ms 
from a 100 Hz low-pass filter shcxild incur the same error as measuring the 
slope at 30 ms from a 50 Hz low-pass filter. 

More important and more subtle is this point: the output of the 

logarithmic electrometer is fairly constant for the initial portion of each 
sweep, then changes to a ramp (slope proportional to l/T^) when electron 
current dominates, and finally saturate in the electron-accelerating portion 
of the sweep. The transition from constant output to ramp occurs fairly 
abruptly for low temperatures (<15 ms for = 400 K) with saturation 
occurring shortly thereafter. For higher temperatures the transition occurs 
much more slowly, the ramp continues for a longer time, and saturation occurs 
later. For example, when * 2000 K the ramp output continues for more than 
90 ms, whereas for = 400 K the ramp output lasts about 20 ms. This 
effect is demonstrated in Figure 6.10. Therefore, lower temper.atures imply 
a shorter period of time between the transition to ramp output and the 
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center of the ramp (where the slope must be measured) . The shorter time 
will cause more error for a given low-pass filter. 

As an exai.iple of error that might occur due to using a low-pass filter 
consider the case T » 400 K with a 450-Hz filter and then again for a 50-Hz 
filter. Suppose the slope is measured from the output signal at t = 10 ms, 
that is, in the center of the 20-ms ramp output. For the 450-Hz filter 
ut » 450 X 2 it X 0.01 = 28. Using equation (6.12) the slope is Y(O.Ol) « a 
and the recovered temperature is 400 K. But for the 50-Hz filter 
tti a 50 X 2tr X 0.01 = 3. The recovered slope now is 0.91 and the apparent 
temperature is 400(1/0.91) = 440 K, an error of 10 percent. Both filters 
give neglible error if t = 45 ms, as would be the case for = 2000 K. 

The above error analysis is in general quantitative agreement with an 
experimental assessment of error due to filtering done in 1967 by L. G. Smith 
(previously unpublished). In the laboratory a 500 K space plasma was 
simulated by using a diode and resistor combination as shown in Figure 6.11. 
The output of the probe circuit was passed through low-pass filters with 
cutoff frequencies between 450 Hz and 20 Hz. The apparent temperature was 
determined, with the error using the 450-Hz filter taken to be zero. The 
results are shown in Table 6.1. Note that the error attributed to using 
a 59-Hz filter is 8.6 percent and that due to a .T5-Hz filter is 11.5 percent, 
for this temperature (500 K) . 

Figure 6.12 is a plot of Y as a function of <ut, as well as ( (a/Y) - IJlOO 
as a function of ut. The latter quantity is the percentage error in 
temperature measurement due to low-pass filtering. 

This concludes the examination of experimental error incurred by 
low-pass filtering in the hardware, both in the paylo.nd and at the telemetry 
station. Errors due to low-pass filtering in the rocket payload appear to 
be ncglibible under all circumstances. It seems that the telemetry low-pass 
filter might conti ibute an error of 10 percent under what may he considered 
a 'worst case' condition (? = 400 K, u) = 50 Hz). The lowest frequency filter 
used on data in this report is 81 Hz; taking Table 6.1 as a guide, it is 
conservative to assume the temperature error in this case is always less 
than 8.6 percent. 

There remain two additional stages of low-pass filtering, these occur 
in computer processing of the data in program ELTF-MP, Within subroutine 


I 



Figure 6.11 Diode and resistor combination used to 
experimentally model 500 K space plasma. 









Figure 6.12 Apparent slope and measurement error as a function of due 
to low-pass filtering. 
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CRUNCH the signal slope i.s determined by a five-^oint least-squares line fit. 
This technique was implemented because of scatter on the data points. The 
resultant slopes still showed significant amounts of variation; to reduce the 
scatter even further the present program averages five consecutive slopes 
to produce a smoothed slope. Both of these operations are low-pass filtering 
operations. However, no error should occur due to this technique because 
even at 400 K the ramp output from the telemetry filter lasts for 20 ms 
(i.e., 20 points) before saturation occurs, whereas only nine points are 
needed to determine the slope. 

Clearly, if too long a period is used to determine the slope (namely, 
a period longer than the duration of the ramp), the effect wi’.l be to 
artifically decrease the slope and elevate the temperature. The nine-point 
slope determination period is short enough to avert any problem in this 
respect. Furthermore, the routine of taking five-point line fits and then 
averaging five slopes yields a smoothed slope equal to 


slope, 




3a. 


3a^ *■ 


2a. 


“5(r 


2*4 - 


3a, 


•Sa^ -2aj 


(6.13) 


where a^, 02 *... , a^ are equally spaced dat.n points. Tins weights the points 
near the center of the interval more than a regiilar nine-point least-squares 
fit would. 


slopC2 


4ag ♦ Sag ♦ 2a^ ♦ la^ - la^ - 2a, - 3a2 - 4Sj 

60 


(6.14) 


This could be advantageous at temperatures below 400 K. Of course, if a 
ramp continues for 9 ms or longer, slo^c, is eqtial to slope,. 

1 4. 

6.4 Probe Surfnac Contnmifuation 

Rocket and laboratory experiments have shown that Langmuir probes nay 
indicate electron energies much higher than present in the ambient plasma if 
the probe has a cont.aminatcd surface [Wehner and Mediaue, 1952; t’irao arid 
Oyama, 1972]. In laboratory studies carried out by Ssuszc&ewios and Holmee 
[1975] they found the electron temperature may he overestimated by as much 
as 100 percent. 

A phenomena log! cal model for the error due to surface cont.aminat ion has 
been proposed by Hrao and Oyama [1972] and critic.iJly tested by Ilirao and 
Oyama in .J.apan and Ssuszc *cwic 2 and Holmes in the United States. The 
contamination layer is modeled as a 'leaky* capacitor in series with the probe. 
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•s in Figure 6.13. Hirao and Oyama estimate the effective values of the 
contamination capacitance and resistancr (for snail probes in the 
E regim) to be on the order of 1 wF and 400 KO. Current to the probe acts 
to charge the contamination layer, so that the potential applied to the 
probe is not the potential seen by the plasma. The first effect of this 
contact potential difference (c.p.d.) is to displace the probe 1-V 
characteristic by an offset voltage. A secondary effect is present for 
Langmuir probes of the type used here. Namely, since the probe dr.iws 
positive-ion current early in the sweep and later draws considerable electron 
current, the voltage to which the contamination layer is charged (c.p.d.) 
will vary during the period of a single sweep and distort the I-V char.ncter- 
istic. Evidently it is this temporal variation that causes error in electron 
temperature measurement. The polarity of the c.p.d. is always such .as to 
effectively reduce the probe sweep rate. Tlic resultant I-V ch.ar.actcristic is 
that for a sweep rate less than that actually applied to the probe. The 
calibration constant is calculated with the actual .sweep r.itc (either 13.5 or 
10.8 V s’^) which results in an overestimation of the electron temperature. 

The contamination layer suggested by Szuszczewicz and Holmes is nothinK 
more than adsorption of neutral atmospheric gas by the probe surface. Tlicy 
provioe the following illustration; consider an initially clean probe 

* 4 

immersed in a neutral gas with u mean molecular weight of 28 amu at 10 imn Hg 

pressure and a temperature of 300 K. These .are values typical of the E region 

at an altitude of 120 km. With unity sticking probability .and a monol.aycr 

14 

defined by approximately 5 x |0 molecules per square cm, the first monolayer 
of contamination develops in approximately 0.13 ms. The next monolayer forms 
in some tens of seconds (due to the decreased probability of sticking) with an 
equilibrium surface condition resulting in sever.al minutes. This contamination 
is in addition to any terrestrial contamination, such as water vapor or finger 
prints. It is obvious that an atomically clean probe is a short-lived 
phenomena, and that attention must be paid to the surf.ace physics of the probe. 

A convenient method of determining when the probe characteristic is being 
adversely affected by surface contamination is to sweep the probe using both 
positive and negative gradients, i.e., *dV/'it .and ’iVjdt. If the character- 
istic is duplicated, then the effect of contamination is small and the ehetron 



Figure 6.15 (a) Effective circuit model for contamination layer as 

proposed by Hirao and Oyama. (b) Hysteresis in con- 
ventional Langmuir probe I-V characteristic resulting 
from probe surface contamination [SsussaseiJtes and Holmes, 

i. •/ / j , 


IS 

quality 
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temperature is correct. If hysteresis is present as in Figure 6.^3’u the 
characteristic is not usable for temperature determination. 

Hirao and Oyama have suggested a technique for avoiding error due to 

varying c.p.d. They sweep the probe at a rate fast enough that the sweep 

period is considerably less than the time constant of the contamination 

layer, R C =0.4 s. Physically, the sweep occurs so quickly that the 
c c 

contamination layer does not have time for its voltage to change 
significantly. In this sweeping mode the effective impedance of the 
contamination capacitance is much smaller than the plasma sheath impedance 
and can be overlooked. Experimentally, a sweep repetition rate of 10 Hz 
was found to do away with surface contamination error. Properly, the 
sweep rate should be given in V s"^, this being the important factor. 

Szuszczewicz and Holmes have adopted a variation of the above technique 
in their development of the pulsed plasma probe. They keep the probe at a 
constant voltage in order to maintain the c.p.d. at a constant value. The 
Langmuir probe characteristic is generated in a point-by-point fashion by 
applying digital voltage pulses to the probe, as shown in Figure 6.14. 

While the voltage pulse is applied to the probe the resultant current is 
sampled. Two conditions must be met for this technique to be effective; 
the pulse period must be much shorter than the layer time constant, and 
much shorter than the interpulse period. A variation of this technique 
would allow the probe to float (drawing no net current) between voltage 
pulses. 

The instrumentation described in this report uses an analog positive- 
going ramp to voltage sweep the probe. Accordingly, there is no hysteresis 
or other effects to indicate if surface contamination is a significant error 
source. With effective repetition rates of only 2 Hz (for the boom probe) 
and 3.3 Hz (for the nose-tip probe), it seems likely that surface 
contamination and the resultant varying contact potential difference could 
be important. In future instrumentation the probe should be alternately 
swept with positive- and negative-going ramps. 
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7. ELECTRON TEMPERATURE IN THE EQUATORIAL E REGION MEASURED 
BY TW) ROCKET EXPERIMENTS AND BY INCOHERENT SCATTER* 


A Nike Apache rocket, 14.532, instrumented with two independent electron - 
temperature probes was launched at 1526 LST into the daytime equatorial 
ionosphere on 28 May 1975 from the Chilca rocket range (12* 30.3' S, 

76* 47.4' W), Peru, as part of the N/.SA operation Project Antarqui. These 
probes were a differential Langmuir probe described by Smith et at, [1968], 
and an ac rectification probe described by Hirao and Oyama [1970]. Personnel 
at Jicamarca Observatory (11* 57' S, 76* 52' W), as co-participants in 
Project Antarqui, simultaneously measured the electron temperatures above 
150 km by incoherent-scatter radar. The ionosphere was quite, with the 
three-hour geomagnetic planetary index (Kp) being equal to 2-. Figure 7.1 
shows the location of the launch site relative to the city of Lima; Jicamarca 
Observatory is northeast from the Lima metropolitan area. 

The operation of the differential Langmuir probe is presented in Figure 
7.2. The upper curve is the typical I-V characteristic for a voltage-swept 
Langmuir probe in the ionosphere. On board differentiation of the I-V curve 
results in the lower curve which is telemetered to ground by standard FM/FM 
telemetry. The electron temperature is recovered by digitally measuring the 
slope of the lower curve between the floating potential and the plasma 
(space) potential. The reciprocal slope in this interval is proportional to 
the electron temperature. 

The operation of the rectification probe is indicated in Figure 7.3. 

The floating potential is a nonlinear function of the amplitude of a 
superposed ac signal and the electron temperature of the ambient plasma. 

The depression of the floating potential is measured for two levels of ac 
voltage which are successively applied to the probe. It is usual to apply 
an ac voltage of a and 2a, the ratio of the floating potential depressions 
being. 



j 2ae 

ts> ^ 


T a§~ 

J 


(7.1) 


*The content of this chapter was presented at COSPAR, 1977 in Tcl Aviv, 
Lsrael, and subsequently published [Smith et al., 1978], 
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Figure 7.1 The location of the rocket range at Chi lea and the 
Jicamarca Observatory are shown relative to Lima, 
Peru. The initial trajectory of Nike Apache 14. .'^32 
is indicated by the dashed line. 
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Figure 7.2 The upper curve is the I-V characteristic for a Langmuir 
probe in the ionosphere. The lower curve, derived from 
the first by on-board differentiation, is used to cal- 
culate r . 
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Figure 7.3 Curve B is the I-V characteristic of the 
rectification probe with no ac voltage 
applied. Curve A is the displaced curve 
due to partial rectification of an applied 
ac voltage; the effect is to depress the 
floating potential. The temperature may 
be determined by measuring the ratio of 
floating potential depressions for two 
levels of applied ac voltage. 
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which may be solved for T , The frequency (20 kHz) of the superposed signal 

is chosen to be lower than the electron*plasma frequency and higher than the 

ion plasma frequency of the ambient plasma. 

The ascending electron-temperature data from both probes are displayed 

in Figure 7.4. The difference, generally being less than 100 K, is attributed 

to experimental error. There is a generally linear increase of electron 

tenqperatui-e from 300 K at 110 km to 1200**K at 180 km. 

Using these data and the experimentally measured electron concentration. 

Figure 7.5, we have made a study of electron heating and cooling in the 

equatorial ionosphere. We used the 1972 COSPAR International Reference 

-22 -2 -1 

Atmosphere, with = 800 K corresponding to 5^^ ^ = 70.6 x lo W m Hz 

for neutral temperature (also shown in Figure 7.4) and con^osition, and the 

discrete cooling formulas given by Datgamo [1969] and Lane and Dalgaxmo 

[1969]. These formulae account for electron cooling due to excitation of 

the atomic oxygen fine structure states, and vibrational and rotational 

states of N 2 and O 2 , We used an oxygen fine structure cooling rate 

coefficient equal to 60% of that suggested by Dalgamo; this is in better 

agreement with a recent cross section calculation by Hoegy [1976]. 

Curves 1 and 2 in Figure 7.6 are the derived cooling rates using T data 

from the rectification probe and the differential Langmuir probe, respectively. 

Curves 3 and 4 are heating rate curves from a numerical model by Dalgaxmo 

et at. [1968] for solar zenith angles x “ 72.1* and 50.2*, respectively. The 

solar zenith angle was 60* at the launch time of Nike Apache 14.532. Above 

140 km the electron heating and cooling rates are in agreement considering 

the uncertainties in the model. Below 130 km the rocket data show cooling 

rates considerably greater than expected on the basis of solar heating, 

suggesting another source of heating. 

As an estimate of the Joule heating at 110 km that might result from 

equatorial electro] et, let us take the current density in the east-west 

-2 

direction to be 10 A km as suggested by Davia et al. [1967] and the 

-3 -1 

conductivity in the east -west direction to be a = 6 x 10 ms m from the 
model by Sugiuna and Cain [1966]. The Joule heating by the electrojet at 
110 km is thus of the order of 1 x 10 ^ eV cm'^ s'\ in agreement with the 
values shown in Figure 7.6. Heating of this magnitude may account for the 
additional source of heating at low altitudes. However, Sahutz et al. [1975] 
using identical instrumentation have observed a similar effect at Wallops 
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Figure 7.5 Electron-concentration profile in the daytime 
equatorial ionosphere. 
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Figure 7.6 A comparison of the calculated electron heating and 

cooling rates. Curve 1 is the cooling rate calculated 
from electron temperature data measured with a recti- 
fication probe. Curve 2 is the cooling rate from 
data taken with a differential Langmuir probe. Curve 
3 is the heating rate from a numerical model by Dalgamo 
for X * 72.1* and Tm * 734 K. Curve 4 is the same 
heating model with x ■ 50,2* and » 766 K. 
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They suggest an additional source of heat below ISO ka due to a 
local aerodynaaic effect. It is clear that if the Joule heating by the 
equatorial electrojet is to be aeasured it is necessary to better understand 
the aerodynaaic effects of a supersonic probe and to have a siaultaneous- 
aeasureaent of the neutral tei^perature. 

During the flight of Nike Apache 14.552, Jicaaarca Observatory obtained 
electron tesqperature data down to 150 ka. These data are displayed in 
Pigi're 7.7 together with the electron te^>erature BMasured by the two rocket* 
borne probes. In the region of overlap (150 to 180 kn) the data are in good 
agreeaent. 


I 



Figure 7-7 Th. average of the electron tenipcratures from 
t :wo rocket -borne probes is shown here 
wiin electron teuperature data obtained 
simultaneously by incoherent scatter at 
Jicamarca Observatory. 
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8. OBSERVATIONS AND DISCUSSIONS 

Electron >t«iq»eraturt profile* for Nike Apache 14.532, 14.534, and 14.533 
are presented in Figures 8.1 through 8.4. Hiere are two profiles for 14.534, 
*>igures 8.2 and 8.3, since data were available using two low-pass filters, 

81 H2 and 450 Hz, respectively. Iliese profiles are plotted on the same axes 
in Figure 8.5 for comparison. 

(hi the aftemoMi of 14.532 there was a constant tm^ierature gradient of 
11*K between 100 and 183 km. These data have been used as the basis for 
a heating/cooling study of the equatorial E region as explained i.". chapter 7. 

(hi the night of 14.534 tha temperature was nearly constant during the 
ascent at about 1350 K or 1300 tC, depending on whether Figure 8.2 or 8.3 is 
taken as correct. This temperature difference is due strictly to the use of 
different low-pass filters in the telemetry station. The 81-Hz filter causes 
the recovered temperature to be elevated by 50 to IOC degrees over the 450-Hz 
data. 

The descent data are between 200 and 300 degrees cooler in both case... 
This is due to a combination of rocket wake effects and, to a lesser degree, 
probe surface contaminants (probe aging). That the wake effect is the major 
error contribution is verified by the tasperature drop occurring rather 
abruptly at apogee and early in the descent leg. A more complete verification 
of this would be to fJy a probe which is not sensitive to probe surface 
effects, such as that described -.y 'Jolntt and SmeaoBooioz 11975). 

The ascent data are considereu valid and representative of conditions 
during a moderately disturbed night. Smith et al. [1974] have considered 
incoherent-scatter data from EVana [1973) and developed the following 
relation for the electron temperature as a function of altitude .and Kp. 

r^(a) • (a - 120) (1.40 Kp ♦ 2.91) * 335*K (8.1) 

This would yield a temperature of 1035 K for an altitude of 20U km and 
Kp • 4. A value of Kp of 6.5 is required in this formula to c.-mse 
temperatures close to 1300 K as observed. It is roughly 300* warmer than 
would be anticipated due to energetic>particle precipitation associated with 
the magnetic disturbance. This is probably due to lessened efficiency in 
one of the electron cooling processes, which would allow a higher final 
temperature to be attained for a given energy flux to the upper .atmosphere. 
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Figure 8.1 Recovered electron-temperature profile for Nike Apache 
14.532. Points connected by the broken line are five- 
point average values. Generally, there is a linear 
temperature gradient with height of about 11 K/km. 
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Figure 8. 
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2 Recovered electron-temperature profile for Nike Apache 
14.534 using an 81-Hz low-pass filter at the telemetry 
station. The points connected by the broken line are 
five-point average values. 
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Figure 8.5 Recovered electron-temperature profile for Nike Apache 
14.534 using a 450-Hz low-pass filter at the telemetry 
station. Points connected by the broken line arc five- 
point average values. 
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Figure 8.5 Comparison of temperature profiles obtained from Nike Apache 
14.534 using different low-pass filters at the telemetry 
station. Heavy low-pass filtering elevates the recovered 
electron temperatures as examined in section 6.3. 
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The temperature profile from Nike Apache 14.533 is much higher than 
can be explained by normal nighttime particle precipitation. Using equation 
(8.1) again, a temperature of 1150 K would be expected for Kp & 5 at an 
altitude of 200 km. Although tonperature data are not available for 
altitudes above 193 km (or on descent) due to an instrumentation failure, 
it se^s reasonable to assume the temperature at 200 km was in excess of 
3500 K. The only explanation presently able to account for the observed 
high tonperatures is that Nike Apache 14.533 passed beneath an SAR arc. 

Battel le-Northwest Laboratory operates a network of computer-controlled 
spectrophotometers which moniter the 6300 A radiation characteristic of SAR 
arcs. On the night of the launch of Nike Apache 14.533 their Michigan 
observatory was inoperative due to hardware problems, and both the Washington 
and Montana observatories were obscured by cloud cover.* Their observatory 
in New York was not yet in operation. There remains some possibility that 
independent verification will come from another ground-based observatory 
or from satellite photographs (although satellite confirmation would be 
questionable unless a 6300 A filter were used, such as on the read line 
photometer aboard the ISIS 2 spacecraft) . Presently, the strongest evidence 
for the existence of an arc that night are the extremely high temperatures 
recorded by Nike Apache 14.533 and the simultaneous observation of a broad 
depression in the P-region electron concentration (density) by the Wallops 
Island ionospheric sounding station. This was noted as characteristic of 
SAR arcs by King and Roach [1961] using bottomside ionograms and later using 
the topside sounder aboard the Alouette II satellite [Clark et al., 1969; 
Norton and Maroviah, 1969; Norton and Findlay, 1969]. 

Wallops Island has a magnetic i-shell value (at 100 km) of 2.6. SAR 
arcs are generally observed between L - 2 and 4 [Roach and Roach, 1963] . 

The mean L-shell value for arcs during the last solar cycle maximum period 
was L = 2.98 [Marovich and Roach, 1963]. 

The thermal gradient between 158 and 187 km is nearly constant at 
41 K km”^. If this is due to downward heat conduction from the base of an 
SAR arc it cannot continue far above 200 km. For example, from visual 
sightings of Sar arcs it has been determined SAR arcs have a well-defined 
base between 300 and 400 km. [Tohnatsu and Roach, 1962; Clark et al., 1969; 

♦Personal communication with D. Slattcr, Rattlesnake Mountain Observatory, 
Battel le-Northwest Laboratory. 
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Roble et oLt 1971]. If the gradient continued to a height of 300 km the 
temperature within the arc would be about 8200 K, whereas temperatures 
measured above arcs by satellites at approximately 1000 km usually peak 
between 4000 and 5000 K [Norton and Findlay, 1969; Chandra et at., 1971; 

Roble et al., 1971; Nagy et al., 1972]. If we take the electron tenqperature 
in the arc to be 4000 K, the observed gradient could continue for only 10 km 
above 187 km. The gradient would then be constrained to remain close to 
zero by the markedly higher thermal conductivity above 200 km. Some evidence 
to this effect might be seen at the top of the profile in Figure 8.4 although 
not enough to substantiate this point. 
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9. RECOMMENDATIONS FOR FUTURE WORK 

The Langmuir probe as used here has an absolute accuracy of ±50 K 
under optimal experimental conditions. There are many sources of error 
and from internal evidence, such as the difference between data obtained 
on ascent and descent of a particular rocket, the error may be ±100 K 
or more. Some improvements seon desirable. 

One possibility is to change the probe sweep to a sawtooth waveform. 
This may allow error due to probe surface contamination to be detected 
and quantified by observing hysteresis effects. 

If surface effects are found to be compromising experimental accuracy 
two routes may be taken: substitution of a faster sweeping rate as 

suggested by Hirao and Oyama [1972], or use of the pulsed -plasma -probe 
technique developed by Holmes and Szuszozewias [1975] . The latter method 
appears to have three advantages. First, the pulse period can always be 
made shorter than the analog sweep period, making the system immune to 
surface effects over a wider range of ionospheric conditions. Secondly, 
data are generated at a more reasonable rate than with fast sweeping, so 
that data buffering is not necessary prior to transmission. And third, 
the pulsed-plasma probe readily interfaces with digital-data transmission. 

Digital -data transmission should be substituted for the analog FM/FM 
telemetry system presently used. This would eliminate telemetry channel 
noise, power line noise, and filtering error. Digital data should be 
directly recorded at the telemetry station with no conversion of any type. 
For a given data rate digital FM transmission will result in a broader 
channel spectrum than analog FM transmission. It will be advantageous 
and reasonable to change the telemetry transmitter modulation index from 
5 to 1 when digital data are being transmitted. 

We suspect the statistical nature of the present electron-temperature 
data [for the nose-tip probe data, anyway) is due to probe surface 
contamination rather than to actual plasma oscillations or plasma 
instabilities due to the vehicle. This was true in Hirao and Oyama’ s study 
[1972]. This should be remedied by a faster sweeping rate of the pulsed- 
probe technique. If the data are still statistical after implementation of 
either of these techniques it would be desirable to generate more data to 
give a higher level of confidence to the data. Accordingly, it might be 
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desirable to separate the electron-concentration ftmction of the probe from 
the electron-temperature function. 

The importance of positive-ion current in electron-temperature 
determination has still not been determined, or proper corrections made. 
Implementation of techniques to do away with surface effects (as recommended 
above) should, at the same time, permit an assessment of the importance 
of positive-ion current. 
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APPENDIX I. Pr«quency response of s logsrithnlc electrometer 

Cmisider the log electrometer circuit shovm in Figure AX.l. A is an 
ideal operational amplifier characterized by infinite input impedance 
(between the * and - terminals), by zero output inqpedance and a voltage gain 
A such that 

V • -Av (I.l) 

The feedback network consists of a diode D, whose properties are 
represented, for forward bias, by 

ijj ■ ig exp{ (u-u J/Wj,} (1.2) 

where ig is the leakage current (for reverse bias) and v^(>ekT/e) is a 
constant with a value of about 25 mV. This formula is valid for 
fp < 10’^ A, by experiment. 

The frequency response of the circuit is limited by capacitance C in 
parallel with the diode. This combines the capacitance of the diode itseh* 
with stray capacitance present in x practical circuit. 

We assume that a current i * nppUed to the electrometer at t b o, 
with y • 0. Then 

f - ip ♦ (1.3) 

where i^ is the displacement current in the capacitor and is given by 

ic - (1.4) 

Using equation (I.l) to eliminate y_ in equations (1.2) and (I 4), 
equation (1.3) becomes 




The solution, for y * 0 at t ■ 0, is 



(1.5) 


( 1 . 6 ) 
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Noting that 4 » 1 and that, aa t t> • v^, aay, the 

equation beeoaes 

V ■ -Wy InjexpC-u^/Uj,) ♦ (1 - exp(-u^/i>pi «*P ^ 

This equation can be put in the paruetric forn 

g-y . g‘Vo ♦ (1 . e^yo) (1.8) 

fdiere y ■ v/v^ and x » (i/^ 2 >)t. The equation is shovm in Figure AI.2 for 

■ 10 and w ■ 20. 

0 

We note that, for x • y 

o 

g~y , g~yo (2 . z Ze'^o or y • y^ - \n 2 (1.9) 

It can be seen in the figure the curves are sinply displaced along the line 
y • X, in agreeaent with eqviation (1.9). 

Again, fron Figure AI.2, the response time may conveniently be specified 
as the tine t at which x u y^, that is 

T ■ Cv^/i (I. 10) 

Froa equation (I.l) )v| » |y_|, since i4 » 1 and 

- y^ ln(i/i^) (I. 11) 

so that equation (I. 10) can be written 
ln(i/ij 

T-CV S- (1.12) 

^ i 

Thus the response tine varies essentially inversely as i. Stated another 
way the frequency response is proportional to i. 

The increase in frequency response as the input current increases is 
illustrated in Figure AI.3. Here the scaled output voltage y ■ v/v^ is 
plotted as a function of scaled time xiji • (i^Cv^ for eight values of 
the input current, expressed by ■ ^^^ 5 * ^^ 1 * 2 er input currents cause a 

proportionally faster change in the output voltage, consistant with the 
conclusion that the frequency response is proportional to input current. 




Figure AI.3 Plot of the output voltage (scaled) of the 

logarithMic electrometer against tu.e (scaled) 
for nine values of input current, ine scaled 
time parameter is xigH ■ where 

is the diode reverse saturation current, C is 
the feedback capacitance, is a constant (25 
mV) and t is the time since the current was 
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It is both interesting and Instructive to compare the preceding 
analysis, which considers a step>£unction increase of current from zero to 
some value i, with the quite-different situation of small changes in current 
about the value i. This is accomplished by considering the dynamic resist- 
ance of the diode. 

From equation (1.2), by differentiation, the dynamic resistance of the 
diode is, with luj » |u 1 (since A » 1) 


Xn ' - ■ ^ 
‘■c 


CI.13) 


where is the average current in the diode. The time constants determining 
the frequency response of the log electrometer for the average current 
is thus 



(1.14) 


The upper 3-dB frequency for the frequency response is 

/„ = — (1.15) 

The frequency response of a log electrometer has been investigated by 

— -8 

Klaus and Smith [1978]. For an average current (i^) of 1 10 A they 

found that the response to a superimposed square wave variation of current 

was a rise time (to 90% of the final value) of ISO ys. This corresponds 

to an upper 3-dB frequency (/^) of 2.5 kHz. 

Using equation (1.14), which = 2S mV we find C = 25 pF. This, then 

is the effective value of the capacitance of the feedback circuit which 

consists of two diodes in parallel, opposed, and stray capacitance. 

The experimental measurements of Klaus and Smith [1978] also show that 

V = 0.24 V for t = 1 X 10 A. From equation (I. 11), therefore, 

^ -12 

= 1 X 10 A. Substituting values for C and i^ in equation (1.12) gives 
the response time for a step-function increase of current from zero. 

T = 6.25 X 10’^^ [ln(10^^t)]/7: 


(1.16) 


This can be interpreted as a frequency response 

/ = « 2.S X lO^Vln(lO^^i) 

aTTT 

-8 

This shows, for example, that f - 500 Hz for i = 2 x lo A. 
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appendix II. Coin)lete listing of Program ELTEMP. Subroutine CRUNCH contains a 
200 ms humbucking algorithm as used in the analysis of Nike Apache 
14.534 (described in section 5.2). 


/JOB 

JOBRRZ. 

SIGNONC 355469629) . 
ZIMMO. 

CHARGE, ELEC, PS6899. 

SETTL,20. 

PRINT/RJEsEE. 

/NOSEQ 


FTN,ER,T,A,Ls0. 

LABEL ( TAPE3 , NT , LB=KU , PO=R , VSNsWI534L-F008 , FsL ) . 

ATTACH , IMSLIB ,UOILIB/UNsLIBRARY . 

♦LIBRARY , IMSLIB , dOILIB . 

LGO. 

SKIPTO,X. 

EXIT. 

ENDSKIP,X. 

DAYFILE. 

/EOR 

PROGRAM ELTEMP ( INPUT , OUTPUT , TAPE 1 = INPUT , TAPE2=0UTPUT , TAPES ) 
INTEGER LOG , ARRAY , RECNUM , START , STOP , WORD , OFFSET 
DIMENSION TEMP(400),D(8) 

LOGICAL FLAG 

COMMON LOG ( 3 , 400 ) , RECNUM( 3 ) , RECTIM( 3 ) , OFFSET , CALIB , KMAX , KMIN 

COMMON/ZIM/ ARRAY ( 2008 ) , NBLOCK , ICHAN 

C0MM0N/Q/T(50) ,ALT(50) ,TIME(400) ,AALT(400) ,IJK,NUM 

COMMON/P/TMLNCH 

DATA ( NBLOCKsO ) , ( SPREAD=250 . O ) 

REWIND 3 

• THE CARD DECK SUPPLIED WITH THIS PROGRAM MUST USE THIS FORMAT. 


CARD 1 ROCKET TYPE AND NUMBER. 3A10 
CARD 2 LAUNCH DATE AND TIME INFORMATION. 3A10 
CARD 3 SOLAR ZENITH ANGLE. 2A10 
CARD 4 LAUNCH LOCATION. 3A10 
CARD 5 DIGITAL TAPE NAME. 1A10 
CARD 6 DIGITAL TAPE RACK NUMBER. 1A10 


CARD 7 TAPE CHANNEL TO BE EXAMINED, II (1,2, 3, 4, OR 5) 

CARD 8 RECORD NUMBER TO START DATA REDUCTION, I4 

CARD 9 RECORD NUMBER TO STOP DATA REDUCTION, I4 

CARD 10 CALIBRATION CONSTANT, F6.0 
CARD 11 OFFSET PARAMETER, 13 (POSITIVE NUMBER). 

CARD 12 KMAX, 14 

CARD 13 KMIN, 14 

CARD 14 LAUNCH TIME DDD;HH:MM:SS 

CARD 15 AND ON SHOULD BE TRAJECTORY DATA SUPPLIED EVERY TEN SECONDS 
THROUGHOUT THE FLIGHT. F0RMAT(F5. 1 ,F'< 1 . 3) 


• INPUT THE CARD INFORMATION. 

READ(1,2)A1,A2,A3,B1,B2,B3,C1,C2,D1,D2,D3,E1,F1 
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2 P0RMAT(3A10/3A10/2A10/3A10/U10/U10) 

RBAD(1,4)ICHAN, START, STOP, CALIB,0FP8BT,KMAX,KMIN,LD,LH»LM,LS 
4 F0RMAT<I1/I4/I4/P5 .0/13/14/14/13 , U, 12 , IX , 12 , IX , 12) 

THLNCHsLD*86400 . 4-LH«3600 . 4-LM060 . •fLS* 1 . 

DO 8 Js1,50 
READ(1,6) T(J),ALT(J) 

6 FORMAT(F5.1,F11.3) 

IF(EOFd)) 10,8 
8 CONTINUE 
10 N0M=J-1 

• OUTPUT THE CARD INFORMATION. 

HRITE(2,12) 

12 FORMAT ("0"/"0",6X,"EEEEEE L",5X,"TTTTTTT EEEEEE M“,5X,''M PPPPP 
J",12X,"AER0N0MY"/" •*,6X,"E",7X,«L'»,8X,"T",5X,"E'*,7X,’'MM MM P 
J PV %6X,"E",7X,"L«,8X,"T'*,5X,'»E'',7X,"M M M M P P",9X, "LABOR 
JATORY"/" ",6X, "EEEEEE L",8X,"T",5X," EEEEEE M M M PPPPP"/" ", 
J6X,"E",7X,"L*,8X,"T",5X,"E",7X,"f« M P",l8X, "ROCKET*/" ",6X,"E 
J",7X,«L",8X,"T",5X,"E",7X,"M M f"/" ",6X, "EEEEEE LLLLLL T 
J EEEEEE M M P",17X,"PROGRAM /"O") 

HRITE(2, 14)A1 ,A2,A3,B1 ,B2,B3,C1 ,C2,D1 ,D2,D3,E1 ,F1 
14 F0RMAT(" ",10X,3A10/" ",10X,3A10/" ",10X,2A10/" ",10X,3A10/" ",10X 
J, "DIGITAL TAPE NAME: ",1A10/" ",1 OX, "DIGITAL TAPE RACK NUMBER: ", 
J1A10) 

WRITE ( 2 , 1 6 ) ICHAN , START , STOP , CALIB , OFFSET , KMAX , KMIN 
16 FORMAT(" ",10X,"TAPE CHANNEL: ",1I1/" ",10X, "START RECORD: ",114/ 
J" ",10X,"ST0P RECORD: «,1I4/" ",10X, "CALIBRATION CONSTANT: ",F6.0 
J/" ",10X, "OFFSET PARAMETER: ",I3/" ",10X,"KMAX: ",I4/" ",10X,"KMIN 
J: ",I4/" ") 

CALL DATE(W) 

CALL GBYTES(H,D(1),6,6,0,8) 

WRITE(2,18)D(4) ,D(5) ,D(7) ,D(8) ,D(1) ,D(2) 

18 FORMAT(" ",10X, "TODAYS DATE: ",2R1 ,"/",2R1 ,"/",2R1/" ") 

WRITE(2,20) (T(J) ,ALT(J) ,J=1 ,NUM) 

20 FORMAT (" ",1 OX, "TRAJECTORY DATA"/" ", 13X, "TIME" ,5X, "HEIGHT", 

J 70(/" ",12X,F5.1,F11.3) ) 

• START UP ROUTINE. 

IC=START-2 
DO 22 J=1,IC 
22 CALL F0RWRD(3,NBL0CK) 

DO 30 N"1,3 
30 CALL LOAD(N) 

• START UP ROUTINE IS COMPLETE. BEGIN DATA PROCESSING. 

IJK=0 
Mz2 
WORDS 1 

35 CALL TRSHLD(M, FLAG, WORD) 

IF (FLAG) GO TO 40 
IF(RECNUM(M).GE.STOP) GO TO 50 
CALL LOAD(LAST(M)) 

MsNEXT(M) 

WORDS 1 
GO TO 35 
40 IJKsIJKi-1 

CALL CRUNCH ( M , WORD , SPREAD , TIME ( UK ) , TEMP ( I JK ) ) 
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HRITE(2,45) TIME(IJK) ,TENP(IJK) 

1»5 FORMAT(“ times", F9. 4," TEMPs",F8.0," DEGREES K."/" ") 

IF(RECNUM(M).GE.STOP) GO TO 50 
CALL LOAD(LAST(M)) 

MsNEXT(M) 

GO TO 35 

• PROBLEM TERMINATION AND DATA OUTPUT. 

50 CONTINUE 

• PARABOLIC ALTITUDE INTERPOLATION 

• RETURNS INTERPOLATED ALTITUDE IN ARRAY AALT(400) 

• FOR BACK TIME IN ARRAY TIME(400) 

CALL TRAJ 

• DATA OUTPUT. 

WRITE(2,55) 

55 F0RMAT("0") 

HRITE(2,12) 

WR1TE(2,14)A1 ,A2,A3,B1 ,B2,B3,C1 ,C2,D1 ,D2,D3,B1 ,F1 
VRITE ( 2 , 1 6 ) I CHAN , START , STOP , CALIB , OFFSET , KMAX , KMIN 
WRITE (2, 80) 

80 F0RMAT("0",T8,«TIME after nominal LAUNCH", T37, 

• "ALTITUDE", T51, "ELECTRON TEMPERATURE"/" ",T15, "(SECONDS)", 
•T35, "( KILOMETERS )",T55,"( DEGREES K)"/" ") 

DO 90 Jsl.IJK 

90 WRITE(2,100) TIME(J) ,AALT(J) ,TEMP(J) 

100 FORMAT(" ",T14,F9.4,T37,F7.3,T59,F5.0) 

STOP 5 
END 

SUBROUTINE LOAD(M) 

INTEGER LOG, ARRAY.RECNUM, OFFSET 

COMMON LOG(3,400) ,RECNUM(3) ,RECTIM(3) , OFFSET, CALIB, KMAX, KMIN 
COMMON/ZIM/ARRA Y ( 2008 ) , NBLOCK , ICHAN 

» SUBROUTINE LOAD IS USED TO PLACE ROCKET DATA IN THE 

• CIRCULAR DATA PROCESSING FILE FROM THE TAPE. THE EFFECT 

• OF A CALL TO LOAD(M) IS TO : 

» FILL ARRAY LOG (M, 1-400) WITH LOG ELECTROMETER DATA. 

• SET RECNUM(M) TO THE RECORD NUMBER OF DATA TRANSFERRED. 

» SET RECTIM TO THE NOMINAL TIME OF RECORD (SECONDS). 

CALL TPGET(3 ,ARRAY, 1 ,NWORDS, NBLOCK) 

IF(NWORDS.EQ.2008) GOTO 20 
WRITE ( 2 , 1 0 ) NBLOCK , NWORDS 

10 FORMATC'O", "RECORD ",I3," CONTAINS ",I4," WORDS.") 

STOP 6 

20 RECNUM(M)=NBLOCK 
DO 30 N= 1,400 
I=ICHAN+N"5 
30 LOG(M,N)=ARRAY(I) 

CALL CALTIM( ARRAY (2006) ,ARRAY(2007) ,ARRAY(2008) ,RECTIM(M) ) 

RETURN 

END 

SUBROUTINE CALTIM( T2006 , T2007 , T2008 , TIME ) 


INTEGER DAY , HOUR ,MIN , SEC , FRCSEC , T2006 ,T2007 , T2008 ,1(16) 
COMfON/P/TMLNCH 

SUBROUTINE CALTIM (CALCULATE TIME) CALCULATES TIMES 
BY DECODING THE NASA DIGITAL TIME DATA IN THE LAST 
THREE WORDS OF A TAPE RECORD. THIS IS UNIVERSAL 
COORDINATED TIME IN DAYS, HOURS, MINUTES, SECONDS, 

AND DECIMAL FRACTIONAL SECONDS. THE NOMINAL LAUNCH 
TIME IS INPUT TO THE PROGRAM IN THE DECK OF DATA 
CARDS. THE TIME AFTER NOMINAL LAUNCH [IN SECONDS] 

IS CALCULATED BY TAKING THE DIFFERENCE BETWEEN THE 
TAPE TIME AND THE LAUNCH TIME. 

CONVERT T2006 INTO FRACTIONAL SECONDS. 

CALL BINARY(T2006,I) 

FRCSECs(8»I(l6)+4»I(15)+2»I(l4)+I(13))»1000 

1+(8»I(12)+4»I(11)+2»I(10)+I(9))»100 

2+(8«I(8)+4»I(7)+2»I(6)+I(5))»10 

3+(8«I(4)+4«I(3)+2*I(2)+I(1)) 

CONVERT T2007 INTO MINUTES AND SECONDS. 

CALL BINARY(T2007,I) 

MINs(4n(15)+2»I(l4)+Kl3))»10 

1+(8»I(12)+4»I(n)+2»I(10)+I(9)) 

SECs(4»I(7)+2»I(6)+I(5))»10 

1+(8»I(4)+4»I(3)+2*I(2)+I(1)) 

CONVERT T2008 INTO DAYS AND HOURS 
CALL BINARY(T2008,I) 

DAY=(2«I(16)+I(15))»100 

1+(8»I(14)+4»I(13)+2«I(12)+I(11))«10 

2+(8»I(10)-f4»I(9)+2»I(8)+I(7)) 

H0UR=(2*I(6)+I(5))»10 

1+(8»I(4)+4»I(3)+2«I(2)4-I(1)) 

CALCULATE THE NOMINAL TIME AFTER LAUNCH. 

TMTAPEsDAY«86400.+HOUR»3600.+MIN»60.+SEC»1 .+FRCSEC* .0001 

TIMEsTMTAPE-TMLNCH 

RETURN 

END 

SUBROUTINE BINARY (IINTGR, I) 

INTEGER 1(16) 

INTGRsIINTGR 

SUBROUTINE BINARY GENERATES THE BINARY REPRESENTATION 
OF THE NUMBER INTGR, FOR USE IN SUBROUTINE CALTIM. 

DO 1 Ns1,l6 

1 I(N)=0 

DO 2 M=1,16 

Js17-N 

L=2«(J-1) 

IF(INTGR.GE.L) I(J)=1 

IF( INTGR. GE.L) INTGR=INTGR-L 

2 CONTINUE 
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• THE ARRAY I(1-l6) IS THE BINARY REPRESENTATION OF 

• THE NUMBER INTQR, WITH I(l6) BEING THE HOST SIQNI> 

• FICANT BIT,2«*15, AND 1(1) THE LEAST SIGNIFICANT 

• BIT,2«0. 

RETURN 

END 

SUBROUTINE TRSHLD(M, FLAG, WORD) 

SUBROUTINE TRSHLDC THRESHOLD) EXAMINES THE LOG ELECTROMETER 
DATA TO LOCATE THE POSITION WHERE SUBROUTINE CRUNCH SHOULD 
BE APPLIED. THE PARAMETERS ARE: 

—M, INPUT. INDICATES WHICH FILE TO EXAMINE. M MAY EQUAL 1,2 OR 3. 
—FLAG, LOGICAL OUTPUT. FLAG IS SET TO .TRUE. IF THE THRESHOLD 
IS FOUND IN THE LOG DATA IN FILE M. 

—WORD, INTEGER INPUT AND OUTPUT. ON INPUT WORD IS A POINTER 
INDICATING WHICH WORD IN FILE M TO BEGIN THE THRESHOLD 
SEARCH. ON OUTPUT WORD IS THE LOCATION OF A BONAFIDE THRES- 
HOLD (IN WHICH CASE FLAG IS SET TO .TRUE.). 

R. K. ZIMMERMAN, JR. MARCH 1978 

AERONOMY LABORATORY - UNIVERSITY OF ILLINOIS 

COMMON LOG ( 3 , 400 ) , RECNUM( 3 ) , RECTIM( 3 ) , OFFSET , CALIB , KMAX , KMIN 
INTEGER LOG, RECNUM, OFFSET, WORD 
LOGICAL FLAG 
FLAGS. FALSE. 

DO 60 JsW0RD,400 

MIsM 

JJ=J+5 

IF(JJ.GT.400) MIsNEXT(M) 

IF(JJ.GT.400) JJsJJ-400 

IF((LOG(M,J).LT. KMAX). OR. (LOG(MI,JJ).GT. KMIN)) GOTO 60 

DO 55 IADDs6,l4 

MIsM 

JJ=J+IADD 

IF(JJ.GT.400) MIsNEXT(M) 

IF(JJ.GT.400) JJsJJ-400 
IF(LOG(MI,JJ).GT.KMIN) GOTO 60 
55 CONTINUE 

• THRESHOLD HAS BEEN FOUND. 

GOTO 70 
60 CONTINUE 

• THRESHOLD WAS NOT FOUND IN FILE M. RETURN. 

RETURN 

• THRESHOLD FOUND • RETURN TO CALLING PROGRAM. 

70 FLAGS. TRUE. 

WORDsJ 

RETURN 

END 

SUBROUTINE CRUNCH ( M , WORD , SPREAD , TIME , TEMP ) 

» SUBROUTINE CRUNCH OPERATES ON THE LOG ELECTROMETER 

• DATA AT THE POSITION INDICATED BY THE POINTER WORD 

• FROM SUBROUTINE TRSHLD. CRUNCH EXAMINES 60 CONSECUTIVE 
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• POINTS, PITTING A LINE TO THE POINTS BY A LEAST SQUARES 

• PROCEDURE (TAKING 5 CONSECUTIVE POINTS AT A TIME). 

• R. K. ZDMERMAN, JR. MARCH 1978 

• AERONOMY LABORATORY - UNIVERSITY OF ILLINOIS 

COMMON LOGO, 400) ,RECNUM(3) ,RECTIM(3) , OFFSET, CALIB,KMAX,KMIN 
INTEGER LOG , RECNUM , OFFSET , WORD , LL0G( 60 ) 

DIMENSION TE(52),SL0PEZ(56) 

• FIRST THE ARRAY LL0G(60) MUST BE FILLED WITH LOG 

• ELECTROMETER DATA. 

N sUORD+OFFSET-t- 1 
NNsHORO-t-OFFSET-i>60 
Ls1 

IF(N.GT.400) GO TO 20 
IF(NN.GT.400) GO TO 40 
DO 10 JsN,NN 
JJsJ-200 
MMsM 

IF(JJ.LT.O) MMsLAST(M) 

IF(JJ.LT.O) JJsJJ<i>400 
LLOG(L)sLOG(M,J)>LOG(MM,JJ) 

10 LxL-i-1 

• ARRAY LL0G(60) IS NOW FILLED WITH DATA FROM FILE M. 

• PROCEED TO LEAST SQUARE FIT. 

GO TO 80 
20 N xN -400 
NNxNN-400 
DO 30 JxN,NN 
JJxJ-200 
MMxNEXT(M) 

IF(JJ.LT.O) MMxM 
IFCJJ.lt. 0) JJxJJ+400 
LLOG(L)xLOG(NEXT(M) ,J)-L0G(MM,JJ) 

30 LxL<»-1 

• ARRAY LLOGC60) IS NOW FILLED WITH DATA FROM FILE NEXT(M) . 

• PROCEED TO LEAST SQUARE FIT. 

GO TO 80 

40 DO 50 JxN, 400 
JJxJ-200 

LL0G(L)xL0G(M,J)-L0G(M,JJ) 

50 LxL+1 
NNxNN-400 
DO 60 Jx1,NN 
JJxJ-200 
MMxNEXT(M) 

IF(JJ.LT.O) MMxM 
IF(JJ.LT.O) JJxJJ+400 
LLOG(L)xLOG(NEXT(M) , J)-LOG(MM, JJ) 

60 LxL+1 

• ARRAY LLOG(60) IS NOW FILLED WITH DATA FROM FILE M AND FILE 

• NEXT(M). PROCEED TO LEAST SQUARE FIT. 


eo DO 170 J>1,56 
INDEX-65 
31*0. 

ST-0. 

SY-0. 

ST2-0. 

SYT-0. 

JP-J-t-M 

IFLAO-0 

90 DO 100 I-J.JP 
RIsI 

RLOOsLLOQ(I) 

S1-SU1. 

ST«ST-i>RI 
SYsSY-fRLOQ 
ST2-ST2fRI*RI 
SYT-SYT+RLOQVRI 
100 CONTINUE 

SL0PE*(S1»SYT-SY»ST)/(S1»ST2-ST»ST) 

B«(SY-SL0PB»ST)/S1 

DO 140 IDUM-1,2 

BADMAXsO. 

110 DO 120 I*J,JP 

IF(I.EQ.IFLAG) GO TO 120 
RI-I 

RLOGsLLOG(I) 

THEOsSLOPE*RI-i>B 

DIFFsABS(THEO-RLOG) 

IF(DIFF.GT.BADMAX) INDEX-I 
IF(DIFF.GT.BADMAX) BAI»1AXsDIFF 
120 CONTINUE 

IF(BADMAX.LT. SPREAD) GO TO 160 
IWORDsWORD-tOFFSET-t-INDEX 
IF(IWORD.GT.400) GO TO 122 
IF(IWORD.LT.I) GO TO 124 
WRITE(2, 126) J, INDEX, RECNUH(M) .lUORD 
GO TO 130 
122 MNXTsNEXT(M) 

IWD4sIWORD-400 

WR ITE ( 2 , 1 26 ) J , INDEX , RECNUM( MNXT ) , IWD4 
GO TO 130 
124 MLST=LAST(M) 

IWD4sIWORD>400 

WRITE ( 2 , 1 26 ) J , INDEX , RECNUM( MLST ) , IWD4 
GO TO 130 

126 F0RMAT("-",15X,''PASS NUMBER*", 12/ 

•" ",15X, "DISCARD POINT (", 12,")"/ 

•" ",15X, "RECORD NUMBER*", 14/ 

•" ",15X,"WORDs",I3) 

130 RlsINDEX 

RLOG*LLOG( INDEX) 

S1sS1-1. 

STsST-RI 

SY=SY-RLOG 


ST2«ST2-Rl»m 

STTBSn-RLOO«III 

SL 0 PB«(S 1 •Srr-SY«ST) / (SI •ST2-ST»SI) 
B«(SY'SL0PB*ST)/S1 
IPLAOiINDBX 
140 COKTINUB 

lfRITB(2,150) 

150 F0RMAT(»-",40X,»1W0 BAD POINTS") 

160 CONTINUB 

IP(SL0PB . LB . 1 . ) SLOPE* 1 . 

3L0P)^Z(J)»SL0PB 
170 CONTINUE 

■ BOXCAR AVERAOINQ OVER 5 CONSECUTIVE SLOPES. 

• COMPUTATION OF TEMPERATURES. 

DO 175 J«1,52 
JP«J-»-4 
SUM>0.0 
DO 173 X>J,JP 
173 SUM«SUM<«^LOPEZ(R) 

SL0PBZ(J)«SUM/5.0 

TB(J)sCALIB/SLOPEZ(J) 

IF(TE(J).QT.3600.) TB(J)«3600. 

175 CONTINUB 

CALL PLOT(LLOO, TE, M, RECNUM) 

CALL MINKTE.TEMP.JFLAO) 

CALL ITIME(JFLAG, WORD, M, TIME) 

RETURN 

END 

FUNCTION NEXT(M) 

IF(M.EQ.3) GO TO 10 
NEXTsM4>1 
RETURN 
10 NEXT*1 
RETURN 
END 

FUNCTION LAST(M) 

IF(M.EQ.I) GO TO 10 
LASTsM-1 
RETURN 
10 LASTx3 
RETURN 
END 

SUBROUTINE PLOT(LLOG , TE ,M, RECNUM) 

• THIS SUBROUTINE PLOTS THE CONTENTS OF ARRAY LLOG AND 

• ARRAY TE FOR EACH EXECUTION OF SUBROUTINE CRUNCH. 

DIMENSION TE(52) >NSYM(2) ,NPTS(2) ,X(60) ,Y(60,2) 
INTEGER LLOG ( 60 ) . RECNUM( 3 ) . IMAG4 (5151) 

DIMENSION A(174),AA(16) 

DATA AA/16»" "/ 
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DATA (AA(1),X«1,8}/«L00 ELECTROMETER VALUE M AND TEM?ERATURB (•) 

* VERSUS NORD NUMBER"/ 

DATA (AA(I),X*9,12)/”NORD NUMBER"/ 

DATA (AA(X),X>13|16)/"L00 VALUE/TEMPERATURE"/ 

DATA (A(X),X«l6l,l6<O/0.0 , 60.0 , 0.0 » 3600.0/ 

DATA (A(X),X«165,166)/"4" , "•"/ 

DATA (NCALLiO) 

XF(NCAa .NE. 0} GOTO 15 
DECODE(160,1000,AA) (A(X) ,X«1 , 160) 

NCALLsI 

15 CONTXNUE 

WRXTE(2,3) M.RBCNUM(M) 

3 FORMAT("0"/"0" 

• ,T32, "//////// RECNUM(",X1,")«",X4,". PLOT(+) XS LOG DAT 
•A, RANGE ZERO ////////"/" ",T32, "//////// TO 1000. (DATA ARE SCALE 
•D TO FXT Y-AXXS). PLOT(«) ////////"/" ",T32, "//////// XS ELECTRON 
•TEMPERATURE, UNSCALED, O.TO 3600.0 . ////////") 

XMXNsO.O 

XMAX»60.0 

• SCALXNO OP LOO VALUES. 

DO 10 Nt1,60 

¥(N,1)s(LLOG(N)-LLOG( 1)455)*3.6000 
XP(Y(N,1).GT. 3600.0) Y(N,1)t3600.0 
IF(Y(N,1).LT.O.O) Y(N,1)«0.0 
X(N)sN 
10 CONTXNUE 
DO 20 N«1,60 
Y(N 2)sTE(N) 

XF(¥(N, 2) .OT. 3600.0) Y(N,2)«3600.0 
XF(Y(N,2).LT.O.0) Y(N,2)«0.0 

X(N)sN 

IF(N.GT.52) Y(N,2)«0.0 
20 CONTXNUE 

CALL USPLX ( X , Y , 60 , 2 , 1 , 60 . A , XMAGU , XER ) 

XF(IER.EQ.O) RETURN 
URXTE(2,30) XER 

30 FORMAT(" ", "ERROR CODE XER * ",I3) 

STOP 7 

1000 P0RHAT(160A1) 

END 

SUBROUTINE MXNI(TE,TEMP, JFLAG) 

DIMENSION TE(56) 

DO 10 1*1,56 

IF(TE(X).LT.250.) TE(I)*50000. 

10 CONTINUE 
20 TEMP*999999. 

JFLAGsO 
DO 30 1*2,55 

IF(TE(I).LE.TENP) JFLAG«X 
IF(TE(I).LE.TEMP) TEMP*TE(I) 

30 CONTINUE 

IFdEMP.GE. 49998.) GO TO 40 


DELTA- (TE(JPLAOf 1)t>TB(JFU0- 1 )) /2-TBMP 
1F(DBLTA.0T.300) TB(JPUO) >50000. 

ZP( DELTA. OT. 300) 00 TO 20 
40 WEITB(2,45) JPUO 

45 FORMATC" »,«•• LOIIEST TEMP 13 TE(",I2,")») 

RETURM 

END 

SUPROUTZNE ITIME(JFLAQ,HORD,M,TIME) 

COmON LOO ( 3 . 400 ) , RECNUM( 3 ) , RECTIM( 3 ) , OFFSET , CALIB , KMAX , KHIN 
INTEGER LOG, RECNUM, OFFSET, WORD 

• R.K.ZI»MERHAN,JR, MARCH 1976 

• AERONOMY LABORATORY > UNIVERSITY OF ILLINOIS 

JUORD> WORD * OFFSET ♦ JFLAO * 4 

• IN THE ABOVE LINE: 

• WORD IS THE WORD IN FILE M THAT QUALIFIED AS THRESHOLD. 

• JFLAO IS THE LOCATION IN LLOO(56) THAT QUALIFIED AS HAVING 

• THE BEST TEMPERATURE WHEN EXAMINED BY MINI. 

• OFFSET IS THE NUMBER OF WORDS BETWEEN THE THRESHOLD WORD 

• AND THE FIRST ENTRY IN ARRAY LLOG. 

• — INTEGER IS TO ALLOW FOR THE FIVE POINT LINE FIT IN 

• CRUNCH. THAT IS, IP LL0G(1) IS INDICATED BY MINI AS 

• HAVING THE BEST TEMPERATURE, THE ACTUAL SLOPE HAS COMPUTED 

• FROM DATA IN LL0G(1), LL0G(2) , . . . . ,LL00(5) . THE TIME 

• ASSOCIATED WITH THE BEST TEMPERATURE SHOULD THEREFORE BE 

• THE TIME OF LL0G(3), THE CENTER OF THE FIVE POINT FIT. 

• THE TIME OP LL0G(3) IS ACTUALLY THE TIME OF L0GG(5), SINCE 

• THERE IS A FIVE POINT BOXCAR AVERAGE OF THE SLOPES BEFORE 

• the FIVE POINT LINE FIT. 

IF(JHORD.GT.400) GO TO 20 
FSEC>(401-JW0RD)«.001 

• NOTE: THE STATEMENT ABOVE CMIPUTES THE FRACTIONAL NUMBER OF SECONDS 

• TO BE SUBTRACTED FROM RECTIM(M) TO GENERATE THE PROPER TIME FOR 

• BEST TEMPERATURE. RECALL THAT THE TIME STORED IN RECTIM(M) 

• IS ACTUALLY THE TIME OF THE FIRST DATA ENTRY IN FILE (M+1 ) . THAT 

• IS WHY WE SUBTRACT BELOW. 

TIME*RECTIM(M)-FSEC 

RETURN 

20 FSEC>(JW0RD-401)«.001 
TIMEsRECTIM(M)4-FSEC 
RETURN 
END 

SUBROUTINE TRAJ 

• SUBROUTINE TRAJ COMPUTES INTERPOLATED ALTITUDE VALUE FOR ANY 

• TIME AFTER LAUNCH. ALTITUDE AND TIME VALUES AT 10 SECOND 

• INTERVALS MUST BE PROVIDED VIA ARRAYS T(50) AND ALT (50). THE 

• ARRAY FTIME(400) CONTAINS THE TIMES FOR WHICH ALTITUDES ARE 

• DESIRED. THESE ALTITUDES ARE CALCULATED AND RETURNED VIA 

• ARRAY AALT(UOO). THE NUMBER OF ENTRIES IN FTIME IS GIVEN BY 

• UK. THE NUMBER OP ENTRIES IN ARRAY T IS GIVEN BY NUM. 
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C0MION/Q/T(50) ,ALT(50),PT1HB(1I00) ,AALT<<iOO) , UK, HUM 

XTCALLiO 

DO 30 J>1,XJK 

XNCRsO 

XP(XTCALL.EQ.I) 00 TO 10 
X«3 

XTCALLil 
5 XNCRsI 
XM1«X-1 
XM2«X-2 

10 XF(FTIME(J).LT.T(XH2)) HRXTE(2,12} 

12 PORMAT("-*',«TXME LESS THAN LOWEST TRAJECTORY POINT.") 

LOCATE THE TXME VALUES HHXCH BRACKET PTXME(J). 

XP(PTXME(J).LE.T(X)) 00 TO 15 

X«l4.1 

XP(X.QT.NUM) STOP 
00 TO 5 

XP PRESENT TXME VALUE PTXME(J) XS IN THE SAME 10 SECOND 
XNTERVAL AS THE PREVXOUS ONE, COMPUTE ALTXTUDE WITH OLD 
COEPPXCIENTS. XP NOT, COMPUTE NEW COEFPXCXENTS. 

15 XF(XNCR.EQ.I) 00 TO 25 
20 AALT(J)« A*FTXME(J)«PTXME(J)<i>B«PTXME(J)4>C 
00 TO 30 

25 BRACU(T(X)-T(XM1))«(ALT(XM1)-ALT(XM2)) 
BRAC2»(T(XM1).T(XM2))«(ALT(I)-ALT(XM1)) 

TOP M BR AC 1 *BR AC2 

BRAC1.(T(XM1)-T(IM2))»(T(X)»T(X)-T(XM1)»T(XM1)) 

BRAC2*(T(X)-T(IM1))»(T(XM1)«T(XM1)-T(XM2)«T(IM2)) 

B0TT0M«BRAC2-BRAC1 

A«TOP/BOTTOM 

Bi(ALTaM1)-ALT(IM2))-A»(T(XM1)«T(XM1)-T(IM2)»T(IM2)) 

B«B/(T(IM1)-T(IM2)) 

C« ALT ( IM2 ) -A»T ( IM2 ) »T ( IM2 ) -B«T { IM2 ) 

00 TO 20 
30 CONTINUE 
RETURN 
END 

SUBROUT INE TPGET ( U . ARRAY , M , N , NBLOCK ) 

THIS SUBROUTINE IS USED TO TRANSFER ONE RECORD OF TAPE DATA FROM 
UNIT «U' TO 'ARRAY', BEOINNINO WITH ARRAY ELEMENT M. THE NUMBER 
OF 16-BIT WORDS TRANSFERRED ON EACH CALL IS N (TYPICALLY 10, 45, 
1005, OR 2008). THE NUMBER OP THE RECORD JUST READ IS NBLOCK. 

IMPLICIT INTBGER(A-Z) 

DIMENSION BUFFER(540), ARRAY(I) 

BUFFER IN (U,1)(BUFPER(1),BUPFER(540)) 

NBL0CK«NBL0CK<^1 
IF(UNIT(U)) 10,20,40 
10 CALL LENGTHX(U,I,J) 

I IS THE NUMBER OF 60-BIT WORDS READ. 

J IS THE NUMBER OF BITS IN THE LAST 60-BIT WORD THAT WERE NOT USED. 
N«(60«I-J)/16 
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CAa 0BYTBS(BUFrtK,AimAy(H)»0j6.0,N) 

RBTURN 

20 NBLOCKtNBLOCA-1 
WRITB(2,30) NBLOCK 

30 PORMATC* "/"O","*** END OP PILE ENCOUNTERED. LAST RECORD IS 

j», §§•«•) 

STOP 1 

AO NRZ1E(2,50) NBLOCK 

50 PORMATC* «/*0*,"»»« PARITY ERROR DETECTED IN RECORD ",IA,*. •••") 
STOP 2 
END 


SUBROUTINE PORURDCU, NBLOCK) 

* *tiIS SUBROUTINE IS USED TO PORWARD SPACE ONE RECORD OK TAPE UNIT U. 

• THE NUMBER OP THE RECORD SKIPPED IS NBLOCK. 

IMPLICIT INTEOER(A-Z) 

DIMENSION BUPPER(SAO) 

BUPPER IN (U,1)(BUPPER(1),BUPPER(5A0)) 

NBL0CKsNBL0CKi>1 
IP(UNIT(U)) 10, 20, AO 
10 RETURN 

20 NBL0CK«NBL0CK-1 
URITE(2,30) NBLOCK 

30 PORMATC "/"O'' ,«•••• END OF FILE ENCOUNTERED. LAST RECORD IS ",IA, 
J*. •••*) 

STOP 3 

AO WRITE (2, 50) NBLOCK 

50 PORMATC" "/"0","»” PARITY ERROR DETECTED IN RECORD •••") 

STOP A 
END 

/EOR 

NIKE APACHE 1A.53A 
1977 AUG 10 05:01:00 UT 
X s 126.5A DEGREES 
WALLOPS ISLAND, VIRGINIA 
WI53AL 
F006 
1 

0925 

1020 

323A0. 

12A 

1000 

0100 


222:05 

:01 :00 

30.0 

026 . 397 

AO.O 

0AA.090 

50.0 

060.730 

60.0 

076. AOA 

70.0 

091.113 

80.0 

10A.87A 
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90.0 117.692 

100.0 129.570 

110.0 140.507 

120.0 150.514 

130.0 159.592 

140.0 167.738 

150.0 174.960 

160.0 181.259 

170.0 186.632 

180.0 191.087 

190.0 194.622 

200.0 197.240 

210.0 198.940 

220.0 199.722 

230.0 199.588 

240.0 198.533 

250.0 196.562 

260.0 193.678 

270.0 189.869 

280.0 185.142 

290.0 179.494 

300.0 172.929 

310.0 165.432 

320.0 157.013 

330.0 147.664 

340.0 137.382 

350.0 126.171 

360.0 114.020 

370.0 100.927 

380.0 086.893 

390.0 071.924 

400.0 056.136 

410.0 039.885 

420.0 023.673 

430.0 008.544 

440.0 0.0 


i 
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APPENDIX III. Program NASA for examining the contents of a digital magnetic 
tape. Sample output for Nike Apache 14.534 is shovm. 


/•CYBER 

JOBROB. 

S1GN0N(355469629) . 

ZIMfO. 

BILL, ELEC, PS6899. 

PRINT/RJEsEE. 

FTN,A,LsO. 

LABEL ,TAPE1 ,NT ,LBsRU ,P0sR , VSNsWI534-D404 , FsL . 

ATTACH , UOILIB/UNsLIBRARY . 

LDSET(LIB=S¥MLIB/UOILIB) . 

LGO. 

SKIPTO,X. 

EXIT. 

OAYPILE. 

/EOR 

PROGRAM NASA (INPUT, OUTPUT, TAPE 1 ,TAPE2sOUTPUT,TAPE3=INPUT) 

• THIS PROGRAM IS USED TO DISPLAY DIGITAL DATA FROM NASA DATA TAPES 

• RECORDED AT WALLOPS ISLAND VIRGINIA. THE USER MUST SUPPLY CARDS 

• INDICATING THE CHANNEL AND RECORD NUMBER TO BE PRINTED. THESE 

• CARDS MUST BE TYPED IN THE FOLLOWING FORMAT: 

• COLUMNS 1 - 4 RECORD NUMBER. 

• COLUMN 5 COMMA. 

» COLUMN 6 CHANNEL TO BE PRINTED - 1,2, 3, 4, OR 5. 

» COLUMN 7 COMMA. 

• COLUMN 8 NUMBER OF CONSECUTIVE RECORDS TO PRINT. 

• FOLLOWING THE LAST CARD IS A 789 CARD. THIS PROGRAM WILL PRINT DATA 

• FOR RECORDS OF LENGTH 10,45,1005, OR 2008 WORDS. EACH WORD IS 16 BITS. 

• WRITTEN BY R.K. ZIMMERMAN, JR. MARCH 1978 

• AERONOMY LABORATORY - UNIVERSITY OF ILLINOIS 

IMPLICIT INTEGER (A-Z) 

DIMENSION ARRAY(2008),PARRAY(400) 

DATA (NBL0CK=0) 

• THE FOLLOWING LINE LOCATES THE 'FILE ENVIRONMENT TABLE' ADDRESS WHICH 

• IS USED LATER IN THE PROGRAM TO BACKSPACE TAPE1. 

NN=FETADR( 1 ) 

• WRITE "NASA" BANNER AT TOP OF OUTPUT. 

WRITE(2,5) 


5 F0RMAT("0"/"0 

",6X 

,"N 


N AAAA 


SSSS AAAA", 13X, "AERONOMY 

J/" ",6X,"NN 

N 

A 

A 

S S 

A 

A" 

J/" ",6X,"N N 

N 

A 

A 

S 

A 

A", 1 0X, "LABORATORY" 

J/" ",6X,"N N 

N 

AAAAAA 

SSSS 

AAAAAA" 

J/" ",6X,"N 

N N 

A 

A 

S 

A 

A", 14X, "ROCKET" 

J/" ",6X,"N 

NN 

A 

A 

S S 

A 

A" 

J/" ",6X,"N 

N 

A 

A 

SSSS 

A 

A", 13X, "PROGRAM" 

J/"0",11X,"— 

CDC 

CYBER 

VERSION — 

"/"O 

«/"0") 


• READ CARD INFORMATION. 

10 READ (3, 20) RECNUM , CHAN , NUM 
20 F0RMAT(I4,1X,I1,1X,I1) 

IF(E0F(3)) 200,25 

• SKIP THE PROPER NUMBER OF RECORDS TO GET TO THE RECORDS DESIRED. 

25 ISKIP=RECNUM-NBLOCK- 1 
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IF(ISKIP) 30,70,50 
30 ISKIPs-ISKIP 
DO 40 Is1,ISKlP 
40 CAa BACKUP(NN,NBLOCK) 

GOTO 70 

50 DO 60 ls1,ISKIP 
60 CALL F0RWRD(1,NBL0CK) 

• FILL ARRAY WITH DATA. 

70 DO 190 Is 1, HUM 

CALL 7PGBT( 1 , ARRAY , 1 ,NWORDS ,NBLOCK) 

IF(NWORDS.EQ.2008) GOTO 90 
IF(NW0RDS.EQ.1005) GOTO 130 
IF(NWORDS.LE.45) GOTO 160 
WRITE (2, 80) NBLOCK,NWORDS 

80 FORMATC* «/"0", "RECORD ",I4," NON-STANDARD LENGTH, ",I4, 

J" WORDS.") 

GOTO 190 

90 DO 100 Ns1,400 
JJsCHAN+N»5 

100 PARRAY(N)sARRAY(JJ) 

CALL CALTIM( ARRAY ( 2006 ), ARRAY (2007) , ARRAY ( 2008 ), DAY, HOUR, MIN, SEC, 
JFRCSEC) 

DlsDAY/100 

D2s(DAY-D1»100)/10 

D3=(DAY-D1«100-D2»10) 

HIsHOUR.'IO 

H2sHOUR-H1»10 

M1sMIN/10 

M2sMIN-M1»10 

S1sSEC/10 

S2sSEC-S1»10 

FlsFRCSEC/1000 

F2s ( FRCSEC-F 1 • 1 000 ) / 1 00 

F3= ( FRCSEC-F 1 • 1 000-F2* 1 00 ) / 1 0 

F4s ( FRCSEC-F 1 • 1 000-F2* 1 00-F3* 1 0 ) 

WRITE (2,110) NBLOCK , ARRAY ( 2006 ) , ARRAY ( 2007 ) , ARRAY ( 2008 ) , D 1 , D2 , D3 , 
JH1 ,H2,M1 ,M2,S1 ,S2,F1 ,F2,F3,F4,CHAN 
110 FORMAK" "/"O", "RECORD ",I4," CONTAINS 2008 WORDS."/" ","TIME CODE 
J (",Z4,";",Z4,";",Z4,") s DAY ",3H,", ".ZII ,":",2I1 ,":",2I1 
J4I1," U.T."/" "."DIGITAL DATA FROM CHANNEL ",I1,".") 

WRITE ( 2 , 1 20 ) ( P ARR AY ( N ) , Ns 1 , 400 ) 

120 FORMAT("0",4X,10I6) 

GOTO 190 

130 WRITE (2, 140) NBLOCK, CHAN 

140 FORMAT (" "/"O" ."RECORD ",I4," CONTAINS 1005 WORDS."/ 

J" "."DIGITAL DATA FROM CHANNEL ", II, ".'•*) 

DO 150 N=1,200 
JJsCHAN+N«5 

150 PARRAY(N)=ARRAY(JJ) 

WRITE ( 2 , 1 20 ) ( PARRAY ( N ) , Ns 1 , 200 ) 

GOTO 190 

160 WRITE(2,170) NBLOCK, NWORDS 

170 FORMAT(" "/"O", "RECORD ",I4," CONTAINS ",I2," WORDS.") 

» THE FOLLOWING THREE DO-LOOPS ARE NECESSARY TO CONVERT TEXTUAL DATA 
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• FROM IBM EBCDIC REPRESENTATION TO CYBER DISPLAY REPRESENTATION. 

DO 173 Ns1,NUORDS 
Ms2*N.1 

173 CALL GBYTES(ARRAY(N),PARRAY(M),44,8,0,2) 

! M«2*NW0RDS 

DO 176 N.1,M 

I 176 PARRAY(N)sCONETD(PARRAY(N)) 

1 DO 178 Ns1,M 

) IF( (PARRAY(N) .EQ.OO) .OR. (PARRAY(N) .EQ.42) .OR. (PARRAY(N) .EQ.50) .OR. 

I J(PARRAY(N).EQ.63)) PARRAY(N)s45 

I 178 CONTINUE 

! HRITE(2,l80)(PARRAY(N),Ns1,M) 

180 FORMATC ",9081) 

190 CONTINUE 

GOTO 10 
200 STOP 
END 

SUBROUTINE CALTIM(T2006,T2007 ,T2008, DAY, HOUR, MIN, SEC, FRCSEC) 
IMPLICIT INTEGER(A-Z) 

DIMENSION I(16) 

• SUBROUTINE CALTIM (CALCULATE TIME) CALCULATES TIMES 

• BY DECODING THE NASA DIGITAL TIME DATA IN THE LAST 

• THREE WORDS OF A TAPE RECORD. THIS IS UNIVERSAL 

• COORDINATED TIME IN DAYS, HOURS, MINUTES, SECONDS, 

• AND DECIMAL FRACTIONAL SECONDS. 

• CONVERT T2006 INTO FRACTIONAL SECONDS. 

CALL BINARY (T2006, I) 

FRCSECs(8»I( 16)+4»I( 15)+2»I( 14)+I(13) )*1000 
1+(8»I(12)+4«I(11)+2«I(10)+I(9))»100 
2+(8»I(8)+4«I(7)+2»I(6)+I(5))*10 
3+(8»I(4)+4«I(3)+2»I(2)+I(1)) 

• CONVERT T2007 INTO MINUTES AND SECONDS. 

CALL BINARY (T2007, I) 

MIN=(4»I(15)+2»I(14)+I(13))»10 

U(8»I(12)+4»I(11)+2»I(10)+I(9)) 

SEC=(4»I(7)+2»I(6)+I(5))*10 

1+(8«I(4)+4»I(3)+2»I(2)+I(1)) 

• CONVERT T2008 INTO DAYS AND HOURS 

CALL BINARY(T2008,I) 

DAYs(2»I(l6)+I(15))»100 

1+(8»I(14)+4«I(13)+2»I(12)+I(11))»10 

2+(8»I(10)+4»I(9)+2*I(8)+I(7)) 

H0UR=(2»I(6)+I(5))*10 

U(8»I(4)+4»I(3)+2«I(2)+I(1)) 

RETURN 

END 

SUBROUTINE BINARY(IINTGR,I) 

INTEGER 1(16) 

• SUBROUTINE BINARY GENERATES THE BINARY REPRESENTATION 


144 


• OF THE NUMBER INTQR, FOR USE IN SUBROUTINE CALTIM. 

INTGRsIINTGR 
DO 1 Ns1,l6 

1 I(N)zO 

DO 2 Ns1,l6 

J=17-N 

L=2«»(J-1) 

IF(INTGR.GE.L) I(J)=1 
IF(INTGR.GE.L) INTGRsINTGR-L 

2 CONTINUE 

• THE ARRAY 1(1-1 6) IS THE BINARY REPRESENTATION OF 

• THE NUMBER INTGR, WITH 1(16) BEING THE MOST SIGNI- 

• FICANT BIT,2»»15, AND 1(1) THE LEAST SIGNIFICANT 

• BIT,2»»0. 

RETURN 

END 

SUBROUTINE TPGET ( U , ARRAY , M , N , NBLOCK ) 

» THIS SUBROUTINE IS USED TO TRANSFER ONE RECORD OF TAPE DATA FROM 

• UNIT 'U* TO ’ARRAY', BEGINNING WITH ARRAY ELEMENT M. THE NUMBER 

• OF 16-BIT WORDS TRANSFERRED ON EACH CALL IS N (TYPICALLY 10, 45, 

» 1005, OR 2008). THE NUMBER OF THE RECORD JUST READ IS NBLOCK. 

IMPLICIT INTEGER (A-Z) 

DIMENSION BUFFER(540), ARRAY(I) 

BUFFER IN (U, 1 ) (BUFFER (1),BUFFER( 540)) 

NBL0CK=NBL0CK+1 
IF(UNIT(U)) 10,20,40 
10 CALL LENGTHX(U,I,J) 

• I IS THE NUMBER OF 60-BIT WORDS READ. 

• J IS THE NUMBER OF BITS IN THE LAST 60-BIT WORD THAT WERE NOT USED. 

N=(60»I-J)/16 

CALL GBYTES ( BUFFER , ARRAY ( M) , 0 , 1 6 , 0 , N ) 

RETURN 

20 NBL0CK=NBL0CK-1 
WRITE (2, 30) NBLOCK 

30 FORMATC ’’/•’O’’,’’*** END OF FILE ENCOUNTERED. LAST RECORD IS ",I4, 
J". •••") 

STOP 1 

40 WRITE (2, 50) NBLOCK 

50 FORMATC "/"O",’’*** PARITY ERROR DETECTED IN RECORD ’’,14,". •••’’) 
STOP 2 
END 


SUBROUTINE FORWRD(U, NBLOCK) 

» THIS SUBROUTINE IS USED TO FORWARD SPACE ONE RECORD ON TAPE UNIT U. 
» THE NUMBER OF THE RECORD SKIPPED IS NBLOCK. 
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IMPLICIT INTEQBR(A-Z) 

DIMENSION BUFFER(540) 

BUFFER IN (U,1)(BUFFER(1),BUFFER(540)) 

NBLOCK-NBLOCK+I 
IF(UNIT(U)) 10,20,40 

10 RETURN 

20 NBLOCKsNBLOCK-1 
WRITE (2 30 ) NBLOCK 

30 FORMAT(" "/"O","*** END OP FILE ENCOUNTERED. LAST RECORD IS •',14, 
J". •••••) 

STOP 3 

40 WRITE (2, 50) NBLOCK 

50 FORMATC ••/"O","*** PARITY ERROR DETECTED IN RECORD ",I4,". •••••) 
RETURN 
END 


INTEGER FUNCTION FETADR(UNIT) 

• THIS INTEGER FUNCTION RETURNS AS ITS VALUE THE ADDRESS OF 

• the "FILE ENVIRONMENT TABLE" CORRESPONDING TO THE FORTRAN 

• UNIT NUMBER GIVEN BY "UNIT", OR ELSE 0, IF THAT UNIT WAS NOT 

• DEFINED IN THE "PROGRAM" STATEMENT IN THE MAIN PROGRAM. FOR 

• INSTANCE, IF THE PROGRAM STATEMENT IS 

• PROGRAM XYZ ( TAPE 1,TAPE2, INPUT, TAPES = INPUT) 

• THEN FETADR(I) WILL BE THE MACHINE ADDRESS OF THE "FET" FOR 

• "TAPE1", FETADR(2) WILL BE ADDRESS OF THE FET FOR TAPE2 

• AND FETADR(S) WILL BE THE ADDRESS OF THE FET FOR FILE "INPUT", 

• SINCE UNIT 5 IS EQUATED TO "INPUT" IN THE PROGRAM. FETADR(6), 

• FOR EXAMPLE, WILL BE 0, SINCE TAPES DOES NOT APPEAR IN THE 

• PROGRAM STATEMENT. 

• THIS FUNCTION IS INTENDED TO WORK ONLY UNDER CDC CYBER FORTRAN. 

• IF YOU TRY TO USE IT IN ANY OTHER CIRCUMSTANCE, YOU DESERVE 

• WHAT YOU GET. ALSO, UNLESS YOU KNOW WHAT A FET IS AND WHAT IT'S 

• FOR, DON'T MUCK WITH THE ROUTINE, BECAUSE OTHER ROUTINES MAY 

• DEPEND ON THIS ONE. 

• CODED FEBRUARY 1978 BY STAN KERR AT THE COMPUTING SERVICES OFFICE 

• OF THE UNIVERSITY OF ILLINOIS AT URBANA, FOR BOB ZIMMERMAN OF 

• ELECTRICAL ENGINEERING. 

IMPLICIT INTEGER(A-Z) 

DIMENSION CORE(I) 

IF(UNIT .GE. 10) GOTO 50 

LFN = 4LTAPE + SHIFT ( 1R0+UNIT, 30) 

GOTO 60 

50 LFN = ilLTAPE + SHIFT (1R0+UNIT/1 0,30) + 

+ SHIFT(1R0+M0D(UNIT,10),24) 

60 B = LOCF(C0RE(0)) 

FETADR = 0 


DO 80 I s 2,64 

IF(C0RE(I-6) .EQ. 0) GOTO 100 
IF((C0RE(I-B) .AND. -777777B) .EQ. LFN) GOTO 90 
80 CONTINUE 

90 FITADR a CORE(I-B) .AND. 777777B 

FETADR a CORE(FITADR-B+l) .AND. 777777B 

100 RETURN 

END 

SUBROUTINE BACKUP ( FETADR , NBLOCK ) 

• THIS ROUTINE IS INTENDED MAINLY FOR BACKING UP 1 BLOCK ON 

• A MAGNETIC TAPE BEING READ BY A FORTRAN PROGRAM. IT IS ASSUMED 

• THAT THE TAPE IS ACCESSED BY THE PROGRAM VIA SOME UNIT NUMBER, 

• SAY 5. THE USER OF THIS ROUTINE SHOULD FIRST USE THE INTEGER 
» FUNCTION "FETADR" TO ESTABLISH THE ADDRESS OF THE "FILE 

• ENVIRONMENT TABLE" FOR THE TAPE, AND SAVE THE VALUE (NEVER MIND 

• WHY, JUST DO IT). THEN THIS ROUTINE MAY BE CALLED AT WILL TO 

• BACKSPACE THE TAPE BY ONE BLOCK. THE SECOND PARAMETER IS SUPPOSED 
« TO BE AN INTEGER VARIABLE WHICH THE USER PRESUMABLY USES TO 

» KEEP TRACK OF WHICH BLOCK THE TAPE IS POSITIONED ON; ALL BACKUP 
» DOES WITH "BLOCK" IS SUBTRACT ONE FROM IT. 

• IF THE TAPE IS UNIT 5, THEN THE PROCESS OF USING BACKUP WOULD BE 

• SOMETHING LIKE THIS; 

• lADDR a FETADR (5) 


CALL BACKUP ( I ADDR , NBLOCK ) 


DO 25 I = 1,10 

25 CALL BACKUP ( lADDR , NBLOCK ) 

ETCETERA , ETCETERA , ETCETERA 

THIS ROUTINE IS INTENDED TO RUN ONLY UNDER CYBER FORTRAN AT THE 
UNIVERSITY OF ILLINOIS. IT USES AN EXTERNAL ROUTINE "SKIPS" FROM 
SYSTEM LIBRARY SYMLIB TO DO THE ACTUAL BACKSPACE. 

CODED BY STAN KERR AT THE COMPUTING SERVICES OFFICE OF THE UNIVERSITY 
OF ILLINOIS AT URBANA FOR BOB ZIMMERMAN OF ELECTRICAL ENGINEERING. 

IMPLICIT INTEGER (A-Z) 

DIMENSION CORE( 1 ) 

BIAS a LOCF(CORE(0)) 

IN = CORE (FETADR+1 -BIAS) .AND. 777777B 

C0RE(FETADR+2-BIAS) a (CORE(FETADR+2-BIAS) .AND. MASK (42)) ♦ IN 


C0BE(FETADRi>3>BIAS) ■ (C0RE(PETADR4>3-BIAS) .AMD. HASK(42)) * IN 


CALL SK1PB(C0RE(PETADR-BIAS) ,1,1) 

NBLOCK ■ NBLOCK-1 

RETURN 

END 


IDENT CONETD 
(( ENTRY CONETD 

•« CONETD - TABLE ORGANIZATION 

• THE CONVERSION TABLE USED BY CONETD IS SET UP WITH 

• 8 DISPLAY CODE CHARACTERS PER WORD, LEFT-JUSTIFIED , TO 

• SIMPLIFIY THE ARITHMETIC NECESSARY TO ACCESS A WORD 

• WHICH CONTAINS THE DISPLAY CODE CORRESPONDING TO A GIVEN 

• EBCDIC CHARACTER. 


CONETD BQ 0 
SA1 XI 


XI s I 


SX6 -1 X6s-1, DEFAULT VALUE IF I OUT OF RANGE 

NG XI, CONETD RETURN IF KO 

SX2 256 

1X2 X1-X2 


PL 

BX2 

AX2 

SA2 

MXO 

BX1 

SX1 

LX1 

BX3 

LX3 

1X1 

SB2 

LX2 

MXO 

BX6 

EQ 


X2, CONETD 

XI 

3 

X2-fTABLE 


57 

-XO»X1 

XU1 

1 

XI 

1 

XUX3 

XI 

B2 

54 

-X0»X2 

CONETD 


RETURN IF I>255 
X2sI/8 

X2 s TABLE WORD WITH CONVERTED VALUE 

X1=MOD(I,8) 

X1sMOD(I,8)+1 


X1s6»(1+MOD(I,8)) 

SHIFT CONVERTED VALUE TO LOW 6 BITS 
AND OUT EXTRANEOUS AND RETURN 


TABLE 


DATA 55626460555255650000B , 55555545564657500000B 
DATA 3334353655555 15500008,4344555561 54737 10000B 
DATA 55555555554742770000B , 55555555550067700000B 
DATA 55554 1 55555555530000B , 55555555374055000000B 
DATA 55555555555555550000B.55556157725145660000B 
DATA 67555555555555550000B , 55556253475277760000B 
DATA 46505555555555550000B.55557556636573710000B 
DATA 55555555555555550000B , 55740060747054640000B 
DATA 55010203040506070000B, 101 15555555555550000B 
DATA 5512131415161 7200000B , 2 1 225555555555550000B 
DATA 55762324252627300000B , 3 1 325555555555550000B 
DATA 55555555555555550000B , 55555555555555550000B 
DATA 72010203040506070000B, 101 15555555555550000B 
DATA 66121314151617200000B.21225555555555550000B 


DATA 7555232M252627300000B , 3 1 325555555555550000B 
DATA 3334353637404 1 420000B , 43445555555555550000B 
BMD 

/EOR 


JOB ACTIVE. 
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— CDC CYBER VERSION — 


RECORD 1 CONTAINS 45 WORDS. 
AMQ1 NIKE APACHE 14.534 


RECORD CONTAINS 10 WORDS. 
AMQ1 


RECORD 3 CONTAINS 45 WORDS. 

AMQ1 NIKE APACHE 14.534 CALS 


RECORD 4 CONTAINS 10 WORDS. 
AMQ1 6 H 


RECORD 5 CONTAINS 1005 WORDS. 
DIGITAL DATA FROM CHANNEL 1. 
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APPENDIX IV. Program CALIB for providing the telemetry calibratim for digital 
magnetic tapes. Sample output for Nike Apache 14.534 is shown. 

/•CtBER 

JOBZIM. 

S1GN0N( 355469629). 

ZIHHO. 

CHARGE, ELEC, PS6899. 

PRINT/RJEsEE. 

FTN,ER.T,A,L80. 

LABEL ,TAPE1 ,NT ,LB»KU ,PO«R , VSN«HI534>0404 , F«L . 

ATTACH , UOILIB/UHsLIBRARY . 

LDSET(LIBsUOILIB) . 

LGO. 

SKIPTO,X 

EXIT. 

ENDSKIP,X. 

DAYFILE. 

/EOR 

PROGRAM CALIB{ INPUT, OUTPUT, TAPE1 ,TAPE2«INPUT,TAPE3*0UTPUT) 

IMPLICIT INTEGER (C*Z) 

DIMENSION CARRAY(2008) ,PARRAY(90) .0(8) 

DATA (NBLOCK«0) , ( INTsO) , (TNI»0) 

REWIND 1 
WRITE(3,10) 

10 FORMAK «0'’/«0“ ,7X , "CCCC" ,4X , ••AAAA'* , 3X , "L" ,7X , "III BBBBB" , 1 2X , 
J«AERONOMY»'/« ",6X,"C C A A L’*,8X,"I B B'*/7X,"C« ,7X, 

J"A A L",8X,'*I B B»,9X,«LAB0RAT0RY«/'* «,6X,'‘C'*,7X,«AAAAAA‘' 

J,» L",8X,»I BBBBB"/" ",6X,"C",7X,"A A L",8X,"I B B", 
J13X,"ROCKET"/7X,"C",4X,"C A A L",8X,»I B B"/" ",7X, 

J"CCCC A A LLLLLL III BBBBB", 13X,"PR0GRAM«/"0") 

CALL DATE(A) 

CALL GBYTES(A,D(1),6,6,0,8) 

HRITE(3,20) D(4),D(5),D(7),D(8),D(1),D(2) 

20 FORMAK" ", "TODAYS DATE: ",2R1 ,"/",2R1 ,"/",2R1) 

READ(2,30) AA,AB,AC 
30 F0RMAT(3A10) 

WR1TE(3,40) AA,AB,AC 

40 FORMAK" ", "TELEMETRY CALIBRATION FOR ",3A10) 

READ(2,50)AD,AE 
50 F0RMAT(A10/A10) 

WRITE(3,60)AD,AE 

60 FORMAK" ", "DIGITAL TAPE NAME * ",A6,", RACK NUMBER s ",A4,«.") 

70 CALL TPGET(1,CARRAY,1,NW0RDS,NBL0CK) 

IF(NWORDS.LE.0045) GOTO 90 
IF(NWORDS,EQ.1005) GOTO 140 
IF(NHORDS.BQ.2O08) GOTO 230 
WR1TE(3,80) NBL0CK,NW0RDS 

80 F0RMAT(" ", "RECORD ",I2," NON-STANDARD LENGTH. ",I4," WORDS.") 

ST0P1 

90 TNIsTNIfI 

00 100 Ns1,NW0RDS 
Ms2*N-1 

100 CALL GBYTES(CARRAKN),PARRAY(M) ,44,8,0,2) 
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M«2«infOI(D3 
DO no N« 1 ,N 

110 PARlUY<lO«COIiSTD(FAmY(N)) 

IP(TNX.BQ.I) UR1TB(3.120) 

120 FORiUT(” *) 

00 125 M«1,M 

IP( (PARRAY(N) .BQ.OO) .0R.(PAiU1AY(N) .BQ.il2) .02.(PARRAY(N) .BQ.50) .OR. 
J(PARRAY(N).EQ.63)) PARRAY(N)«A5 
125 CONTZIIUB 

HRZTBC3.130) NBL0CK,(PARRAY(N),Nb1,N) 

130 PORMAT(" "."RBCORD ".I2.1X.90R1) 

GOTO 70 
no INTsZIITfl 

ZP((I1IT.BQ.1).AND.(TNZ.EQ.0)) VRZTS(3,150) AD 
150 P0RMAT(*’0”,”N0 ID HEADER RECORDS ON TAPE ".Ae,”.”) 

SUM1 bSUM2«SUM3bSUM4bSUM 5»0 
DO 160 Hs5, 1000,5 

N1sN+1 $ N2sM-«>2 $ N3 bN 43 $ NMsN^A $ N5«N+5 
SUH1>SUMUCARRAY(N1) 

SUM2«SUM2<»CARRAY(N2) 

SUN3bSUN3-*-CARRAY(N3) 

SUHA«SUMAi>CARRAY(NA) 

160 SUH5«SUM54CARRAZ(N5) 

AltSUMI/200. $ A2«SUM2/200. A A3tSUH3/200. $ A4sSUMA/200. 
A5>SUM5/200. 

IP(INT.EO.I) HR1TE(3,170) NBL0CK.A1 ,A2,A3,A4,A5 
1P(INT.EQ.2) VRITB(3.180) NBL0CK,A1 ,A2,A3,A4,A5 
IP(INT.BQ.3) VRZTE(3,190) NBL0CK,A1 ,A2,A3,A4,A5 
IP(INT.BQ.A) VRITE(3.200) NBL0CK.A1 ,A2,A3.AH,A5 
IP(ZNT.BQ.5) WRZTE(3.210) NBL0CK.A1 ,A2,A3.A4,A5 
IP(ZNT.QE.6) URZTB(3,220) NBL0CK.A1 ,A2,A3,A4,A5 

goto 70 

170 r0RMATO*0'*,26X,«CHAN 1 CHAN 2 CHAN 3 CHAN A CHAN 5**/ 

J " "."RECORD ",Z2," LBE (-7.505)",5F9. D 

LBH (-3.755)",5F9.1) 

BAND CENTER ".5F9.D 
UBH (♦3.75»)",5F9.1) 

UBE (♦7.50*)",5F9.1) 

".5F9.1) 

230 IP((TNX.GE.1).AND.(ZNT.EQ.O)) WRITE(3,240) AO 
1F((TNI.EQ.0).AND.(INT.EQ.0)) WRZTE(3.250) AD 
240 FORMAT ("0", "NO TELEMETRY CALIBRATION RECORDS ON TAPE ",A6,".") 

250 FORMAT ("0", "NO HEADER ID RECORDS OR TELQfETRY"/ 

J " "."CALIBRATION RECORDS ON TAPE ",A6.".") 

URITE(3.260) NBLOCK.AD 

260 P0RMAT("0". "RECORD ".12." IS THE FIRST DATA RECORD ON TAPE ". 
JA6.".") 

ST0P4 

END 


180 FORMATC" "."RECORD 
190 FORMAT(" "."RECORD 
200 FORMAT(" "."RECORD 
210 FORMAT(" 

220 FORMAT (" 


"."RECORD 

"."RECORD 


M2." 

". 12 ." 

",I2." 

",I2." 

",I2." 


SUBROUTINE TPCET ( U . ARRAY . M . N . NBLOCK ) 


• THIS SUBROUTINE IS USED TO TRANSFER ONE RECORD OF TAPE DATA FROM 

■ UNIT 'U* TO ’ARRAY*. BEGINNING WITH ARRAY ELEMENT H. THE NUMBER 


• OP 16-BIT HOBBS TBANSPSBBBO OM BACH CAa IS M (TYPICAaY 10, A5, 

• 1005, OB 2006). THB NUHBBB OP THB BBCOBD JUST BBAD IS MBLOCK. 

IMPLICIT INTBOBB(A-Z) 

DIHBMSION BUPPBB(540), ABBAY(I) 

BUPPBB IN (U,1)(BUPPBB(1),BUPPBB(5B0)) 

NBLOCK-NBLOCX^I 
IP(UNIT(U)) 10,20,60 
10 CALL LBN0THX(U,I,J) 

• I IS THE NUHBBB OP 60-BIT HOBDS BBAO. 

• J IS THE NUHBBB OP BITS IN THB UST 60-BIT HOBD THAT WEBB NOT USED. 

N«(60*I-J)/16 

CALL 0BnBS(BUPPBB,ABBAY(H),0,l6,0,N) 

BETUBN 

20 NBLOCK*NBLOCK-1 
WBITE(2,30) NBLOCE 

30 F0BMAT(" «/'*0'',''»«» END OP PILE BMCOUNTEHBD. LAST BECORD IS ",I6 

jii, •••••) 

STOP 1 

40 HRITE(2 50) NBLOCK 

50 FORMAT(« "/"O","*** PARITY ERROR DETECTED IN RECORD ".IM,". •••" 
STOP 2 
END 


SUBROUTINE FORWRD(U, NBLOCE) 

• THIS SUBROUTINE IS USED TO FORWARD SPACE ONE RECORD ON TAPE UNIT U. 

• THE NUIffiER OF THE RECORD SNIPPED IS NBLOCE. 

IMPLICIT INTEGER (A-Z) 

DIMENSION BUFFER (540) 

BUFFER IN (U,1)(BUFFER(1),BUFFER(540)) 

NBL0CX>NBL0CE4>1 
IF(UNIT(U)) 10,20,40 
10 RETURN 

20 NBLOCKsNBLOCK-1 
WRITE(2,30) NBLOCE 

30 FORMATC" END OF FILE ENCOUNTERED. LAST RECORD IS ",I4 

J". •••") 

STOP 3 

40 WRITE (2, 50) NBLOCK 

50 FORMATC" "/"O","*** PARITY ERROR DETECTED IN RECORD ",I4,". •••" 
RETURN 
END 

IDENT CONETD 
( ( ENTRY CONETD 

•« CONETD - TABLE ORGANIZATION 

• THE CONVERSION TABLE USED BY CONETD IS SET UP WITH 

• 8 DISPLAY CODE CHARACTERS PER WORD, LEFT-JUSTIFIED, TO 

• SIMPLIFIY THE ARITHMETIC NECESSARY TO ACCESS A WORD 


o o o o o 
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f 

CONBTD 


TABLE 


CCCC 

C 


C 

CCCC 


NMiai CONTAIKS THE DZSPUT COOS CORKBSPONDINO TO A OZVBN 
EBCDIC CHARACTER. 

EO 0 

3A1 XI 

SX6 -1 

NO XI, CONBTD 

SX2 256 

ZX2 X1-X2 

PL X2, CONBTD 

BX2 XI 

AX2 3 

SA2 X24>TABLB 

MXO 57 

BX1 -X0*X1 

SX1 XU1 

LX1 1 

BX3 XI 

U3 1 

XXI XUX3 

SB2 XI 

LX2 B2 

MXO 54 

BX6 -X0«X2 

EQ CONETD 

DATA 55626460555255650000B.55555545564657500000B 
DATA 33343536555551550000B,43445555615473710000B 
DATA 55555555554742770000B , 55555555550067700000B 

DATA 55554155555555530000B,55555555374055000000B 
DATA 55555555555555550000B , 55556 1 57725 1 45660000B 

DATA 67555555555555550000B.55556253475277760000B 
DATA 46505555555555550000B , 5555755663657 37 1 0OOOB 

DATA 55555555555555550000B.55740060747054640000B 
DATA 55010203040506070000B, 101 15555555555550000B 
DATA 551213141516 1 7200000B , 2 1 225555555555550000B 

DATA 55762324252627300000B.31325555555555550000B 
DATA 55555555555555550000B , 55555555555555550000B 

DATA 72010203040506070000B, 101 15555555555550000B 
DATA 6612131415161 7200000B , 2 1 225555555555550000B 

DATA 75552324252627300000B.31325555555555550000B 
DATA 33343536374041420000B.43445555555555550000B 
END 
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ROCKET 
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LLLLLL 

III 

BBBBB 

PROGRAM 


XI ■ I 

X6i-1, DEFAULT VALUE IP I OUT OF RANOB 
RETURN IP KO 


RETURN IP I>255 
X2iI/8 

X2 s TABLE HORD WITH CONVERTED VALUE 

X1sH0D(I,6) 

X1«K0D(Z,6)<»>1 


X1«6«(UM0D(I.6)} 

SHIFT CONVERTED VALUE TO LOW 6 BITS 
AND OUT EXTRANEOUS AND RETURN 


TODAYS DATE: 07/19/78 

TELEMETRY CALIBRATION FOR NIKE APACHE 14.534 
digital tape name = HI534 , RACK NUMBER = D404. 


RECORD 1 AMQ1 NIKE APACHE 14.534 

RECORD 2 AMQ1 

RECORD 3 AMQ1 NIKE APACHE 14.534 CALS 


RECORD 

4 

AMQ1 6 H 





CHAN 1 

RECORD 

5 

LBE (-7.50J) 

168.0 

RECORD 

6 

LBH (-3.75») 

1101.0 

RECORD 

7 

BAND CENTER 

2033.8 

RECORD 

8 

UBH (+3.7511) 

2967.2 

RECORD 

9 

UBE (+7.50$) 

3903.3 


CHAi^ ' 

CHAN 3 

CHAN 4 

CHAN 5 

151.0 

169.1 

154.1 

166.9 

1098.8 

1106.8 

1088.3 

1113.3 

2039.6 

2050.8 

2030.0 

2055.1 

2986.6 

2985.8 

2971.8 

2998.: 

3932.0 

3931.3 

3911.2 

3944.1 


RECORD 10 IS THE FIRST DATA RECORD ON TAPE HI534 
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APPENDIX V. Subroutine CRUNCH as used in the analysis of Nike Apache 14.533. 

Seventeen points are averaged to smooth the data and remove 60 Hz hum. 

SUBROUTINE CRUNCH (H, WORD, SPREAD , TIME , TEMP) 

• SUBROUTINE CRUNCH OPERATES ON THE LOG ELECTROMETER 

• DATA AT THE POSITION INDICATED BY THE POINTER WORD 

• FROM SUBROUTINE TRSHLD. CRUNCH EXAMINES 60 CONSECUTIVE 

• POINTS, FITTING A LINE TO THE POINTS BY A LEAST SQUARES 

• PROCEDURE (TAKING 5 CONSECUTIVE PCINT.S AT A TIME). 

• R. K. ZIieiEimAN, JR. MARCH 1978 

• AERONOMY LABORATORY - UNIVERSITY OF ILLINOIS 

COMMON LOG ( 3 , 400 ) , RECNUM( 3 ) , RECTIM( 3 ) , OFFSET , CALIB , KMAX , KMIN 
INTEGER LOG , RECNUM , OFFSET , WORD , LLOG ( 60 ) 

DIMENSION TE(52),SL0PEZ(56) 

• first the ARRAY LLOG (60) MUST BE FILLED WITH LOG 

• ELECTROMETER DATA. 

N =WORD+OFFSET+ 1 
NN=WORD-i-OFFSET-i-60 
L=1 

• SEVENTEEN POINT AVERAGE. 

T17*0.0 

NXsWORD-iOFFSET 
DO 5 NZIM=1,17 
JJsNZIM-9-t-NX 
MMsM 

IF(JJ.LT.I) GOTO 3 
IF(JJ.GT.400) GOTO 4 
T17sT17+LOG(MM,JJ) 

GOTO 5 

3 JJsJJ+400 $ MM=LAST(M) 
t1T=T17+L0G(MM,JJ) 

GOTO 5 

4 JJsJJ-400 $ MM=NEXT(M) 

T17=T17+L0G(MM,JJ) 

5 CONTINUE 

IF(N.GT.400) GO TO 20 
IF(NN.GT.400) GO TO 40 
DO 10 JsN.NN 
JLSS=J-9 $ JPLS=J+8 
MLSS=MPLS=M 

IF(JLSS.LT.I) MLSS=LAST(M) 

IF(JLSS.LT.I) JLSS=JLSS+400 
1F(JPLS.GT.400) MPLS=NEXT(M) 

IF(JPLS.GT.hOO) JPLS=JPLS-400 
T 1 7=T 1 7-LOG ( MLSS , JLSS ) +LOG ( MPLS , JPLS ) 

LLOG(L)=T17/17.0 
10 L=L+1 

• ARRAY LLOG (60) IS NOW FILLED WITH AVERAGED DATA FROM FILE M. 

• PROCEED TO LEAST SQUARE FIT. 
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GO TO 80 
20 N sN -400 
NNsNN- 400 
DO 30 JsN.NN 
JLSSsJ-9 $ JPLSsJ<f8 
MLSSsMPLSsNBXT(M) 

IF(JLSS.LT.I) MLSSsM 

IF ( JLSS . LT . 1 ) JLSSs JLSSi-400 

T 1 7 a T 1 7 -LOG ( MLSS , JLSS ) +LOG ( MPLS , JPLS ) 

LL0G(L)aT17/17.0 
30 LaL+1 

• ARRAY LLOG(60) IS NOW FILLED WITH AVERAGED DATA FROM FILE NEXT(M) . 

• PROCEED TO LEAST SQUARE FIT. 

GO TO 80 

40 DO 50 J=N,400 

JLSSaJ-9 $ JPLSaJl>8 
MLSSaMPLSaM 

IF(JPLS.GT.400) MPLSaNEXT(M) 

IF ( JPLS. GT. 400) JPLSaJPLS-400 
T17aT17-LOG(MLSS,JLSS)+LOG(MPLS,JPLS) 

LL0G(L)=T17/17.0 
50 LrL+1 
NNaNN-400 
DO 60 Ja1,NN 
JLSSaJ-9 $ JPLSaJ+8 
MLSSsMPLSsNEXT(M) 

IF(JLSS.LT.I) MLSSaM 

IFCJLSS.lt. 1 ) JLSSaJLSS-i-400 

T 1 7 =T1 7-LOG ( MLSS , JLSS ) +LOG ( MPLS , JPLS) 

LL0G(L)aT17/17.0 
60 L=L+1 

• ARRAY LL0GC60) IS NOW FILLED WITH AVERAGED DATA FROM FILE M AND FILE 

» NEXT(M). PROCEED TO LEAST SQUARE FIT. 

80 DO 170 Js1,56 

THE REMAINDER OF SUBROUTINE CRUNCH IS IDENTICAL 
WITH THAT GIVEN IN APPENDIX II. 
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